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I. INTRODUCTION 


This is the fifth and final annual technical report for NASA Grant 
NGR-11-002-133 (March 1971 - April 1976) "Investigation of Lubricant 
Rheology as Applied to Elastohydrodynamic Lubrication". A brief 
summary of the first four years of research is included in Section II. 
This report is primarily concerned with the research performed during 
the last year. 

During the first four years significant progress was made in 
three areas; a) rheology of lubricants in a high pressure rheometer, 

b) observation of lubricant behavior in an elastohydrodynamic simulator, 
including the development of a temperature measuring technique, and 

c) the analytical modeling of elastohydrodynamic contacts (isothermal 
non-linear rheological fluid model and a linear model with thermal 
effects) . 

During the past year the research was concerned with elasto- 
hydrodynamic simulation and exploration of glassy state lubricant 
behavior. The ehd simulator research (Section III) is a continuation 
of the infrared temperature measurements to study the maximum surface 
temperatures in ehd contacts including rolling contacts with sliding. 

The effect of roughness on surface temperature and wear has also been 
studied. Correlations of surface temperature and operating variables 
have been developed which are suitable for use in engineering design. 

Preliminary studies of the glass transition of lubricants were 
undertaken to determine if typical lubricants in typical ehd contacts 
pass through a glass transition (Section IV) . They clearly indicate 
that some lubricants do pass into the glassy state in ehd contacts. 
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These results require further study and will result in a change in 
thinking regarding ehd lubricant rheology. 

As a necessary part of the glassy state studies an analytical pro- 
cedure was developed which relates the measured ehd film thickness to 
the pressure distribution. No assumption about lubricant rheology was 
necessary. This effort is presented in Section V. 

Appendix B presents a list of publications which resulted totally or 
in part from this research. 

Some of the research described in this report was jointly supported 
by NASA and the National Science Foundation (Grant Number ENG 74-21002). 

W. 0. Winer and D. M. Sanborn were the principal investigators and 
had overall responsibility. Others were involved in specific aspects 
of the research: M. A. Alsaad conducted the light-scattering glass 

transition experiments, which constituted his doctoral research; S. Bair 
conducted the other glass transition experiments and assisted in the ehd 
experiments; H. S. Nagaraj conducted the ehd temperature measurement 
experiments, which constituted a major portion of his doctoral research; 
and G. M. Rentzepis and R. K. Kunz were responsible for the elasticity 
analysis of Section V. 
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II. SUMMARY OF RESEARCH CONDUCTED DURING 1971-1975 


A brief summary of research completed during the first four 

years of this grant is presented in this section. For more complete 
details in each area the previous annual contract reports or referenced 
publications should be consulted. 

A. High pressure Rheology of Lubricants 

Fifteen fluids have been subjected to high pressure rheological 
measurements. These include several fluids which consist of pairs of 
base fluid and base fluid plus an antiwear additive. These measurements 
were in two categories: 1) determination of the effect of pressure and 
temperature on the fluid viscosity (and pressure-viscosity coefficients) 
at relatively low shear stress (approximately lkPa [1,2], and 2) determination 
of the effect of high shear stress upon the viscosity at modest pressures 
and temperature [3,4]. The pressures and temperatures chosen for the latter 
experiments were those which gave a low shear stress viscosity comparable 
to the average apparent viscosity observed in the ehd contact simulator 
under the sliding point contact conditions employed in this laboratory. 

All fifteen fluids have been subjected to the former type of measurements 
while six fluids have been subjected to the latter. 

The pressure viscosity measurements have shown that the pressure- 
viscosity coefficient, important in ehd, decreases with increasing tempera- 
ture in such a manner that the natural logarithm of the pressure-viscosity 
coefficient is approximately inversely proportional to the temperature [1,2]. 
Maximum differences of less than a factor of four in pressure-viscosity 
coefficient were observed among typical lubricants at a temperature of 
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37. 8C (100 F) . These differences are reduced considerably at 149C (300 F) . 
Differences in viscosity and viscosity-pressure coefficients between base 
fluids and base fluids formulated with antiwear additives were found to 
be minor in all cases [1,2]. 

The high shear stress measurements performed on six fluids were 
conducted up to a shear stress within a factor of 3 to 4 of the average 
shear stress occurring in the sliding point contact ehd simulator employed 
in this laboratory (5 MPa compared to 14-20 MPa) . Repeatable data 
explainable in terms of viscous heating only were obtained on four of the 
six fluids throughout the entire range of pressure, temperature and shear 
stress examined. The other two fluids, a polymer blended mineral oil and 
a 50 cs dimethyl siloxane, exhibited rheological behavior which could not 
be explained by viscous heating alone. The polymer blended mineral oil 
exhibited shear thinning behavior at relatively low shear stress as expected. 
Both the polymer blended mineral oil and the dimethyl siloxane exhibited 
"non-liquid" behavior above some characteristic threshold pressure and 
shear stress. In both cases the apparent viscosity during the non-liquid 
behavior was higher than expected from Newtonian considerations and 
generally was neither repeatable nor correlatable with the chronology of the 
measurement. In retrospect this behavior may be related to glass transition 
in these materials at the conditions of the experiments. 

During the high shear stress measurements the average resident time 

of a fluid particle in the capillary, hence the time it was subjected to 

the high shear stress, was short (4 ps) compared to the average time the 

2 

fluid particle is in an ehd contact (10 ps) . Although the capillary 
residence time varied from a minimum of 4 ps up to 10^ ps no effect of 
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residence time on viscosity was observed for the four fluids mentioned 

above suggesting time dependent effects are not important in this range. 

This is consistent with the oscillatory crystal research of Lamb and 

co-workers which indicates that the shear relaxation time for fluids 

-3 -2 

in this viscosity range are of the order of 10 to 10 ys. 

Viscous dissipation is an important rheological phenomena in the 
high shear behavior of typical liquid lubricants. This phenomena can be 
confused with shear thinning pseudoplastic behavior. The utilization 
of capillaries of different length-to~diameter ratios in performing the 
experimental measurements is essential in distinguishing between pseudo- 
plastic and viscous heating effects. In addition to interpreting the 
experimental data, an analysis of the flow of a Newtonian fluid in a 

capillary including viscous dissipation and temperature dependent viscosity 
was performed [5]. The results of this analytical model are consistent with 
the experimental measurements. The analytical model was able to predict 
the high shear stress behavior of the lubricants whose behavior was measured 
with a good accuracy up to the highest shear stress measured. 

Jakobsen's film thermal-shearing model [6] has been applied to point 
contacts [7] with emphasis on comparing predictions with existing data from 
this laboratory. The model appears to have some merit in conditions of high 
sliding. However, one result of this study has been the recognition of the 
possible importance of lubricant behavior in the glassy state. This 
possibility is discussed in Section V of this report. 

B. Elastohydrodynamic Simulator Studies 

All the fluids that were studied in the high pressure rheometer have 
also been studied in the ehd simulator where film thickness distribution 
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and traction were measured. Most of these measurements were at Hertz 
pressures in the range 1.0 - 2.0 GPa and with sliding speeds in the range 
0.02 - 5.0 m/s. 

The most significant advancement in these studies has been the 
development of the infrared emission technique for determining the fluid 
and ball surface temperature distributions in the ehd contact [8,9,10]. 

The technique utilizes an indium-antimonide liquid nitrogen cooled 
detector to measure the IR radiation emitted by the contact region. By 
suitable filtering, the radiation from the ball surface and the fluid 
can be separated and with the system calibration each temperature can be 
determined. At lGPa Hertz pressure the temperature measurement area 
resolution is 1% of the contact area. Two principal applications of the 
technique have been studied. The first was to measure the temperature 
distribution for a given fluid and then employ these temperatures in an 
analytical model assuming a given fluid constitutive equation to predict 
traction which was compared with the measured value [ 7 ]. Because of the 
rheometer measurements mentioned above the linear constitutive model with 
temperature and pressure dependent viscosity was employed. The second 
major application of the IR temperature measuring technique was in observing 
the behavior of an ehd film as the severity of the contact condition 
increases [9,10]. The objective was to learn about failure of severely loaded 
ehd contacts. Based on the surface temperature mappings that have been 
completed to date, temperatures in some regions of the contact are higher 
than expected when considering viscous dissipation as the only energy 
source. The analytical model discussed below permits the prediction of 
the maximum fluid temperature through the film for a given film thickness. 
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Based on this model the maximum predicted temperature would be in the 
center of the contact. However, the measured maximum temperature may be 
in the lobe region and considerably higher than predicted for some 
operating conditions [9], The region where measured temperatures exceed 
the predicted temperatures coincided with the region where the optical 
film thickness was equal to or less than the composite surface roughness 
of the surfaces utilized [9], Thus both the film thickness and temperature 
measurements suggest the possibility of asperity interaction resulting 
in energy dissipation in addition to that due to viscous dissipation. 

Throughout the experiments performed thus far, ball surface tempera- 
tures have been obtained. In addition, the fluid film temperature and ball 
surface profile have been determined at selected points. The ball surface 
temperature variations with speed are in close agreement with the results 
predicted by Jaeger, Blok, and later Archard, for temperature variations 
of run in surfaces [10], These correlations are further developed in 
Section III of this report. Ball temperatures as high as 250 C for a reser- 
voir temperature of 35 C have been observed. The lubricant temperatures 
have been approximately 40 C above those of the ball surface. 

C. Analy tical Modeling of Ehd Contacts 

Two analytical models of ehd contact behavior have been developed. 

One was for the complete solution of the isothermal lubrication of rollers 
with a non-linear viscous fluid covering the load range from the lightly 
loaded rigid cylinder case to the highly loaded elastic cylinder case [11], 
The second model is concerned only with the prediction of traction in highly 
loaded elastic contacts with viscous dissipation and temperature dependent 
linear viscosity [6,7]. 
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The complete solution for the isothermal lubrication of rollers 
with non-linear fluids has been a difficult numerical problem to execute. 
However, solutions are now available with both linear and non-linear 
fluids for lightly loaded rigid cylinders and for elastic cylinders [11,12], 
in which case the solutions have been completed to the limit of the computiqg 
capabilities available. When this research was initiated it was thought 
that shear thinning of lubricants at high shear stress may be an important 
rheological phenomena. However, the concurrent work mentioned above of 
the high shear stress measurements of viscosity now indicate that shear 
thinning probably is not important for most lubricants. Those findings 
make the results of this analytical modeling less important than originally 
expected. The results are, however, still expected to be of value for two 
reasons; a) the numerical technique employed is different from most 
previous ehd analytical studies and therefore the results can serve as a 
comparison with other works for the linear fluid case, and b) because of 
the broad similarities between shear thinning and the effect of viscous 
heating the results of this model may give some qualitative insight into 
ehd behavior where viscous heating is present. The results will of 
course be useful for applications employing non-linear shear thinning 
fluids, however, this class of lubricants is now believed to be a 
relatively small portion of liquid lubricants currently employed. There 
may be wider application of the results to those systems employing 
greases, however. 

The second analytical model developed [6,7] was the result of the 
high shear stress rheometer measurements which indicated the importance 
of viscous heating compared to shear thinning at high shear stress. This 
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model is concerned with predicting traction in the contact given a 
velocity, film thickness distribution and a rheological constitutive 
model for the fluid. The constitutive model employed is a linear model 
with pressure and temperature dependent viscosity. The model also 
accounts for the thermal effects of viscous energy dissipation. Because 
of the dependence upon obtaining film thickness independent of the analysis 
this effort is primarily useful in assessing the validity of rheological 
constitutive relations. Preliminary application of the model to a sliding 
ehd contact suggests that a linear viscous model with temperature dependent 
viscosity and viscous dissipation may be sufficient to describe the 
behavior of the fluid in an ehd contact with high sliding. However, the 
model fails to predict the traction as the amount of sliding decreases. 

This failure of the model may be the result of the lubricant behavior as a 
glassy solid in low sliding and rolling ehd contacts. 
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III. EXPERIMENTAL INVESTIGATIONS IN EHD CONTACTS 
A. Introduction 

The surface temperature data obtained for the case of pure 

sliding under the previous contract [ 13 ] was analyzed and compared 

to predicted flash temperatures. The surface temperature data is 

in good agreement with the predicted flash temperature. 

Both the surface roughness profile and the AC signal from 

the infrared detector were analyzed in the frequency domain. It was 

found that only certain wavelengths in the roughness profile are 

important as far as asperity interaction is concerned. The parameter 
* 

A(A = h c /a) is shown to be a valuable means of predicting surface 
interaction in EHD contacts. A new method for determining the surface 
roughness, cr, is also proposed. 

The effect of various amounts of rolling on the film thickness, 
traction and surface temperature was studied and the results are 
reported. A correlation study of the influence of sliding and rolling 
velocity on surface temperature was also undertaken and the results 
are discussed. 

•flp 

The nomenclature used in this section is found on the following 

pages . 
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NOMENCLATURE 


a 


C 

d 

h 

i 

h 

i 

k 


L 

£ 

P(u>) 

P H 

q 


r 



t 


T 

AT 


T 

o 

TC 


V 

V 


V 


BP 


V 

s 


W 


Hertzian contact radius, m 

Auto-correlation function n 

Hertzian contact diameter = 2a ,m 
Film thickness at contact center, m 
Minimum film thickness , m 
Thermal conductivity, W/ (miK) 

Non-dimensional time parameter = 

Length, m 
Power spectrum 

2 

Peak Hertz contact pressure, N/m , Pa 
2 

Heat flux, W/m 
Correlation coefficient 
Coefficient of determination 
Time, s 

Thermal diffusion time = a^a,s 

Characteristic resident time = a/V,s 

Temperature, C 

Temperature rise, C 

Upstream (bulk) temperature, C 

Traction coefficient 

Surface Velocity, m/s 

Rolling velocity = 1/2 (V^ + V ga > , m/s 

Break-point velocity in correlation study, m/s 

Sliding velocity = - V ga , m/s 


Load , N 
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X 


Distance along surface, m 

2 

a Thermal diffusivity , m /s 

3 Length shift in computing autocorrelation function, m 

* 

3 Correlation distance, m 

K A constant = .968 a 1/2 a 1/2 /k, m 3/2 s 1/2 C/N 

A h /a 

c 

A Wavelength, ra 

2 

JJ Viscosity, Ns/m 

T. Slide roll ratio = V /V 

s 

O Surface roughness, rms or Ra, m 

T Time , s 

Wave numbers = 1/A, m ^ 

Subscripts 
b ball 

f fluid 

sa sapphire 
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B. Experimental Equipmen t and Techniques 


Two EHD contact simulators were used in the studies conducted 
during the past year. The first (Figure III-l) is the apparatus used for 
film thickness, traction and temperature measurements of EHD contacts 
in pure sliding. It is the same as the equipment used in pure sliding 
EHD studies previously conducted in this Laboratory. The new 
equipment (Figure III-2) differs from the sliding rig primarily in that 
the sapphire disk may be rotated by a separate drive system to achieve 
various amounts of slip. The new equipment also has a simplified 
system for traction measurement. In the pure sliding apparatus, the 
traction force was isolated from the normal force by using an air 
bearing under the sphere support. In the new rig, the lateral motion 
needed to detect the traction force is obtained by mounting the ball 
support on sheet metal columns. The system is therefore very rigid 
with respect to the normal load, but is relatively flexible with respect 
to the traction force. Strain gages are attached to the columns such 
that bending strain (proportional to traction) can be detected, but 
compression (proportional to normal load) can not. 

In each simulator, the contact is formed using a 31.8 mm diameter 
52100 steel ball rotating and loaded against a sapphire flat. Lubricant 
film thickness was measured for the smooth ball experiments using an 
optical interference technique. Both monochromatic and dichromatic [14] 
methods were employed. 

Chromium steel balls (AISI 52100) of three different roughness 
values have been used to study the effects of surface roughness. The 
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Ra values are 0.011, 0.075 and 0.38 ym and are referred to as smooth, 
medium rough and rough, respectively. These balls do not have a 
preferred roughness orientation. Because of the high loads and relatively 
rough surfaces used in this study, it was not always possible to accurately 
measure film thickness by an optical technique. The fluid temperature 
can not be obtained without an accurate film thickness distribution. 

Rough surfaces give a time dependent film thickness at a given location 
which prohibits the determination of film temperatures for the rough 
ball. Therefore, only ball surface temperatures were measured. 

The technique for measuring the ball surface temperature 
distribution in the contact has been described earlier [9,12], Ball 
surface temperature measurements are made onlv at the contact center 
for the case of smooth and medium rough balls. The center of the EHD 
contact is easy to locate and also, the temperature at this location 
is quite close to the maximum. In the experiments with the rough 
ball, the temperature measured is the maximum in the contact, whether 
or not it occurs at the contact center. It was necessary to make this 
procedural change because the contact center could not be located when 
the rough ball was used. This is due to the disappearance of the inter- 
ference fringe pattern because of light scattering at the relatively 
large surface asperities. Therefore, the contact was scanned for the 
maximum surface temperature in each case. 

The infrared radiation emitted at this contact is measured with 
an infrared radiometric detector having a spot size resolution of 38 ym 
with a 15X objective and a response time of 8 ys in the A.C. mode of 
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operation. It was necessary to monitor the oil reservoir temperature with 
a thermocouple and to hold it constant with a constant temperature bath and 
an oil circulation system. 

A Bendix surface profilimeter has been equipped with a relocation 
device [13]. This allows the same area of the ball surface to be 
examined in the profilimeter both before and after use in the EHD simulator. 

Fourier analysis of both the surface profile data and the high 
frequency I-R data was done using Hewlett-Packard Model 5451A Fourier 
Analyzer System. This equipment is owned by another unit of the Institute 
but has been made available for use on his project. 

The fluid used throughout this investigation is a naphthenic 
base oil designated N1 in previous studies. A complete description 
of the fluid is given in the appendix of reference [15]. 

C . Flash Tempe rature Measurements for Pure Sliding 

Ball surface temperatures for the smooth ball were reported a 

year ago [13], but a correlation study with existing flash temperature 

theories has only just been completed. Figure III-3 shows a plot of the 

ball temperature rise at the contact center above the bath temperature, 

as a function of sliding velocity for peak Hertz pressures ranging 

2 

from 0.52 to 2.03 GN/m . For a given Hertz pressure, the data 
plotted on log-log coordinates falls on one straight line for velocities 
up to a break-point value and then on another line, of lower slope, 
for higher velocities. The temperature rise is the difference between 
the ball surface temperature at the contact center and the lubricant 
bath temperature as measured with a thermocouple. 

The trends shown in Figure III-3 can be predicted using the 
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techniques of Blok [16], Jaeger [17] and Archard [18]. The energy 
dissipation rate in the contact is equal to the product (W»TC-V^). 

The heat flux to each surface is then proportional to the dissipation 
rate divided by the .area of the contact, or 


W-TC-V, 

q * ^ “ TC *V V b 


( 1 ) 


"ira 

where the constant of proportionality includes the portion of the 

total heat flux transferred to the ball. According to Archard [18] 

the temperature rise on the surface of a solid moving at velocity 

and subjected to heat flux q can be described in terms of a non- 

2 

dimehsional parameter L = t^/t 2 * The term t^ = a /2 ot^ represents 
the time required for the effect of q to penetrate a distance (a) 
below the surface, whereas t^ = a/V represents the time required for 
a point in the contact to move a distance (a). Therefore, for a given 
fluid and surface, L is proportional to (P *V, ) . The variation in AT 

n D 

(ball surface temperature minus bulk oil temperature) is given as 


AT oc q.p n a TC * V • P 

H b H 


for L < .1 


( 2 ) 


and 


- 1/2 - 1/2 
AT « q-p/ -V 


TC-V^-P 2 ' 2 

D n 


for L > 5 


(3) 


The parameter L has been plotted in Figure III-3 for the case of P„ = 1.70 

rl 

2 

GN/m . The values 0.1 and 5 represent points for which the full solution 
is in good agreement with the limiting cases of equation 2 (stationary 
heat source solution) and equation 3 (fast moving heat source solution) 
respecitvely . Most of the data shown in Figure III-3 has an L between 
the limiting values. 
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In order to determine how well the data fits the Jaeger-Ar chard 
predictions of equations (2) and (3) , a multiple regression analysis 
of the data shown in Figure III-3 was performed. For each Hertz 
pressure level, the data was divided into two regions separated by 
a break-point velocity clearly visible in Figure III-3. This break- 
point velocity is a function of Hertz pressure and can be described by 

V - °- 53 V 1 (4) 


with a correlation coefficient of r = 0.95. This results in a value of 
3.4 for L at the break point in the experimental data. This is reason- 
ably close to unity which is the theoretically predicted breakpoint [18]. 
The reason for this discrepancy is the limited amount of experimental 

data at low values of L. For sliding velocities below V-,,-, 

BP 


AT 


2.02 0.53 

H 


(5) 


with R 


2 


.978. 


For velocities 
AT - P r 


above V Bp , 
1.34 v 0.34 


( 6 ) 


with R 2 = .985. 

These relations were obtained neglecting the traction coefficient, 
since traction data was not obtained during these experiments. However, 
from equations 1-3 it is apparent that its influence must be included. 
Based on previously reported data [19,20], the traction coefficient has 
only a slight positive dependence on P in the pressure ranges studied. 
Furthermore, TC was found to be dependent on sliding velocity to the 
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power -0.29 with r = .994. The data used had 3 £ L £ 22. Including 
this dependence in equations 5 and 6 results in 


AT « p 2 - 02 v* 82 V- 29 -c p 2 * 02 v- 82 
H b b H b 


TC 


(7) 


and 


AT cc p^*^ V * 29 cc p^-"3^ 

H n b H b 


TC 


( 8 ) 


Recalling that equations 2 and 3 represent limiting cases, the experi- 
mental results appear to be in good agreement with the predicted trends 
(equations 2 and 3) . 

D . Correlations Between Surface Roughness and Temperature 

In the previous report [13 ] , a method was described wherein the 
surface roughness profile could be determined at the same location on 
the ball surface both before and after running in the EHD simulator. 
During the past year, a significant extension of this work has been 
undertaken on the frequency or wavelength aspects of the surface profile 
signal and its influence on surface temperature. 

In order to study the ball surface profile in the frequency or 
wavelength domain, the profilometer output was first recorded on 
magnetic tape. This was done for samples of balls from the three 
roughness classifications, both before and after running in the EHD 
simulator . 

The recording of the profilometer trace was transferred to a 
Hewlett-Packard Fourier Analyzer. The first step in any analysis 
performed by this instrument is an A/D conversion. After analysis 
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of the digital representation of the input signal, the results may be 
printed and/or a D/A conversion can be made and the results plotted 
or displayed on a CRT. 

Figure III-4a shows the surface profile of the unused ball 

having an initial surface roughness of 0.38 pm Ra. The trace represents 

the signal after the A/D and D/A conversions and is an excellent 

replica of the original signal plotted by the Bendix profilometer . The 

distance along the ball surface is shown in terms of the Hertzian 

contact diameter, d. For the conditions of steel on sapphire at 1.24 

2 

GN/m peak Hertz pressure, the calculated contact diameter is 0.43 mm. 

Figure III-4b shows the same ball after running for 24 minutes 

in the simulator. The profile traces shown in Figure III-4 were taken 

in the direction of the sliding motion. No attempt was made to measure 

the profile at precisely the same location on the ball surface. The 

traces taken are of sufficient length that they should be representative 

of any position within the wear track. From the previous report [13], 

2 

the conditions of 1.24 GN/m and V = 1.02 m/s should result in 

s 

a value of A < 1 based on the unused ball surface roughness. The value 
A < 1 indicates that asperity interaction should be expected. Comparing 
the two traces in Figure III-4, it is apparent that the 24 minute run-in 
period has resulted in a smoother surface. The roughness is .25 pm Ra 
which results in A - .4. The 15.9 mm radius of curvature of the ball 
surface can also be seen. 

1. Power Spectra 

The primary function of the Fourier Analyzer is to perform the 
Fast Fourier Transform (FFT) on the time varying innut sipnal . A conveni 
ent way of displaying the amplitudes of the component terms in Fourier 


19 



series is through the power spectrum. The input to the analyzer is a 

time varying voltage. After the A/D conversion, the FFT results in a 

complex amplitude for each of the 127 (plus a constant) terms used in 

the series. Each of these terms is assigned a pair of channels in the 

analyzer. The frequency range, and therefore, frequency increment can 

be selected. The frequency assigned to two adjacent pairs of channels 

differs by a constant. The resulting plot of the power spectrum is 

linear in the independent variable, frequency. The value of the 

ordinate at each center frequency (channel) is the product of the complex 

amplitude for that frequency and its conjugate. The units are, therefore, 
2 

V versus Hz. The D/A conversion makes the plotted results appear to be 
continuous, but all of the data processing was done digitally in terms 
of the coefficients of the first 127 terms of the Fourier series. 

Although the results are plotted on a linear frequency scale, 
this is not a convenient representation for this study. This is due 
to the fact that the frequency is influenced by the profilometer stylus 
scanning speed and differences between magnetic tape recording and 
play-back speeds. In order to remove these effects, the abscissa 
has been changed to read wavelength. The unit of wavelength 
selected is d, the Hertz contact diameter of 0.43 mm. The power 
spectra of the profiles shown in Figure III-4 are given in Figure III-5. 

The first four channels (0-3) were cleared prior to plotting. 

This is because the values are so large compared to the values in higher 
channels (see insert) , that a plot would show only the very low 
frequency components. Therefore, only components having wavelengths 
<3d have been plotted. In comparing the power spectra for the unused 
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and run-in ball surface profiles, note the factor of two change of 
scale in Figure III-5a versus III-5b. The most apparent contrast is the 
reduction of the components having wavelengths <2d. Comparing the 
tabular data in Figure III-5, it appears that roughness components having 

wavelengths >2d remain essentially unchanged after the 24 minute run-in 
period. The nine percent difference in DC level can be attributed to 
slight differences in signal amplification in obtaining the two surface 
profiles. The most drastic change occurs for a wavelength of one 
Hertz diameter. This component has been reduced by more than a factor 
of eight. 

2 . Autocorrelation Function 

In addition to the power spectrum, the autocorrelation function 
can be useful in describing the surface profile. In fact, Whitehouse 
and Archard [21] and Peklinik [22] suggest the description of a real 
surface in terms of only G, the rms roughness, and B*, the correlation 
distance. Briefly, the autocorrelation function is obtained by com- 
paring a time varying signal, in this case the surface profile measure- 
ment, with a replica of itself where the replica is shifted some amount 
of time (or distance at a given velocity) . If the shift along the sur- 
face profile is 3, then the autocorrelation function C(3) becomes 


C(B) = -q f a(x)a(x 4- 3)dx 


where a(x) is the amplitude of the surface profile trace at a distinct 
x along the surface. For a truly random signal, C(3) will be a maximum 
at 3 = 0. If the signal is periodic, C(3) will peak whenever 3 is a 
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multiple of the wavelength. If a Gaussian distribution of surface 
roughness amplitudes is considered, equation 9 results in an exponential 
autocorrelation function when normalized such that C(0) = 1.0, 

C(B) = exp (-3/ 3*) (10) 

Therefore, when 3 = 2.3$*, C(B) = 0.10. The value 0.10 has been 
arbitrarily set by Whitehouse and Archard as being sufficiently small 
that two points on the surface may be regarded as being independent. 

The power spectrum is the inverse Fourier transform of the autocorrela- 
tion function, or 


P(w) = Owe ia)T dT (11) 

Substituting equation 10 into 11 results in 

P(w) = 0.3 for A = 1/go = 3* 

P(go) - 0 for A = 1/go < 0.13* (12) 

P(go) - 1 for A = 1/co > 103* 

The wavelength 3* can therefore be thought of as a dividing point 
between those surface wavelength components having significant ampli- 
tudes and the shorter wavelengths which do not. Whitehouse and Archard 
have shown the model to be in good agreement with a large number of real 
surface profiles. 

Figure III-6 shows the autocorrelation function for the unused 
and run-in ball surfaces shown in Figure III-4. The position shift 3 
is shown in terms of the Hertz diameter d. If the Whit'ehouse-Archard 
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model is applied to the data in Figure III-6, the value of 8* can be 
found from C(2.38*) = 0.1 C(0). The correlation distance 8* is ,17d 
(74ym) for the unused ball and .33d (142ym) for the run-in surface. 

It is apparent that the exponential model of equation (10) can only 
be used for 8 £ 2.38*- The further increase in C(8) implies that the 
profile is not truly random, but has strong periodic components with 
wavelengths > d. This can also be seen from the power spectra shown 
in Figure III-5. 

Figure III-6 also shows that the correlation distance has doubled 
as a result of the run-in process. From equation 10, a doubling of 8* 
implies that two points on the surface must be taken twice as far apart 
after run-in to be considered independent. Because of the interrela- 
tionship between the power spectra and autocorrelation function (equa- 
tions 11 and 12) , if the surfaces were in fact random, the half power 
point would be shifted to a wavelength twice the value of that for the 
unused surface. This is the case when the high frequency components 
are reduced much more than the low frequency components (see Figure 
III-4) . This is also consistent with the notion of independent points. 
Consider, for example, a surface comprised of only two wavelength 
components 2d and 0.2d. When the 0.2d wavelength component is present, 
the autocorrelation function would have a value of 1.0 for 8 = 2md 
where m is an integer. However, it would also have sizable peaks 
for 8 = 0.2md (with drops to near zero between). After run-in, 
however, the 0.2d wavelength component will be substantially reduced. 
Assuming it to be eliminated, the C(8) now peaks only once every 
8 = 2md and reaches zero between the peaks. If the correlation distance 
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is defined as the ten percent point, it is apparent that this distance 

must be increased if high frequency components are reduced or eliminated. 

3. Correlation between Surface Rouehness and Temperature 
Fluctuations 

The AC signal from the infrared detector obtained during the 

run-in process was also recorded and analyzed on the Hewlett-Packard 

analyzer. However, in order to eliminate noise in the detector signal, 

a 20 KHz, 3dB filter was used during data acquisition. The noise 

~6 2 

remaining below 20KHz had a peak power of about 5 x 10 V whereas 

-4 2 

the detector signal had a power content on the order of 10 V . 
Furthermore, the sampling technique used by the Fourier analyzer tends 
to average out the random noise while enhancing the signal. 

Figure III-7 shows a normalized histogram of the AC signal from 
the detector obtained at the beginning of the run-in period (Figure 
III-7a) and after running for eight minutes (Figure III-7b) . The figure 
shows that the voltage distributions are nearly Gaussian. Secondly, it 
is apparent that the amount of dispersion is decreased by the run-in 
process. It is interesting that both the high and low extremes in 
temperature signal have been eliminated as the asperity peaks were 
removed. The symmetry of the fluctuations supports the argument 
advanced previously [10] that the peak in the temperature fluctuation 
signal is due to local heating at an asperity; whereas the minimum 
signal occurs directly after this maximum and is due to the relatively 
low local pressure in the region following an asperity. 

Figure III-8 shows the power spectra for the same IR detector 
signals as used in Figure III-7. The abscissa has been replotted in 
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terms of wavelength in units of d (Hertz diameters) . It is apparent 
from Figure III-8 that the run-in process results in a significant 
reduction in power at all wavelengths. In both spectra, the peak 
occurs near A * 2d, although the run-in process has reduced the peak 
magnitude by more than a factor of three. It should be noted that the 
detector signal was also filtered to remove frequency components less 
than 500 Hz. This corresponds to wavelengths in Figure III-8 which 
are ereater than 5d. 

4. Discussion of Results of Spectra Correlations between 
Temperature Fluctuations and Surface Roughness 

Based on the results shown in Figures III-5 and III-8, the 
authors believe that the methods presently used to determine 
acceptable levels of the Ra or rms surface roughness may not be 
the most appropriate when surface roughness effects on EHD 
contact temperature, wear or scuffing are of interest. For example, 
the rms value for the surface shown in Figure III-4a (unused rough 
ball) is equal to the square rool of the integral of the power 
spectrum (Figure III-5a) . However, the proper limits of integration 
must be chosen. The Bendix instrument used to obtain the profile 
shown in Figure III-13a uses an upper limit of .76 mm (through 
electrical filtering) and has a lower limit set by the stylus radius 
of 13ym. If one assumes that the absence of surface temperature 

fluctuations is indicative of negligible asperity interaction (as 
shown in the previous report [13]), then as far as asperity inter- 
action is concerned, one can determine the significant wavelengths 
of the surface profile. For example. Figure III-8a shows that the 
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temperature spectral power content during asperity interaction 
drops to about ten percent of its peak value for the wavelengths 
A d/4. However, Figure III-5a shows that there are significant 
roughness components down to at least A = .08d (the power content is 
essentially the square of the amplitude) . The same figures show that 
the major components for both the roughness profile and the temperature 
signal are in the range d £ A £ 2d. Therefore, whereas the roughness 
wavelength in the range d £ A 1 2d are obviously important in causing 
temperature fluctuations, wavelengths shorter than d/4 appear to have 
no effect. This, of course, is based on the premise that a temperature 

fluctuation is caused by an asperity interaction. Based on the 
power reduction to ten percent of its maximum value, the lower bound 

for rms determination should be A - d/4. The upper limit is set at 
A = 2d because the temperature signal drops off rapidly for A > 2d 
even though there are roughness components present with wavelengths 
greater than 2d (see Figure III-5) . 

On the basis of the data taken, it appears that the rms surface 
roughness to be used in a calculation of A(A = h/d) should consider 
only wavelengths in the range d/4 £ A £ 2d. For this experiment, 
d = 0.43 mm and the applicable range is then .11 1 A £ .86 mm. 

A roughness of about . 38ym rms was obtained using the Bendix system 
with limits of about .013 £ A f .76 mm, and the standard filter 
(2CR-type). However, by integrating the power spectrum of the surface 
profile in the Fourier Analyzer between the limits of 0.11 £ A £ 0.86 mm 
the rms value for the surface roughness was 0.53yra. This higher value 
is believed to be due to significant power in range 0.76 £ A £ 0.86 mm. 
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The Bendix instrument also has available cut-offs at .25 mm and .067 mm. 
The use of either of these values would correspond t ms integration 
limits of . 030d _ _ . 58d and ,030d - _ .18d respectively. Since 

the desired range is .25d - _ 2d , the resulting rras values from the 

profilometer using these cutoffs would not be expected to be a 
reliable measure of surface roughness in predicting surface interaction. 

In addition to the possibility of getting a value of surface 
roughness which cannot reliably predict performance, the fact that 
wavelengths shorter than 0.25d do not appear to be significant might 
have an economic impact. For example, under conditions in which the 
Hertz diameter is relatively large, the final stages in a lapping or 
polishing process may not be important. 

The EHD contact appears to behave like a mechanical filter in 
which wavelengths outside of a relatively narrow band have a negligible 
effect on the performance of this system. If relatively thin films 
are expected during the operation of a highly loaded bearing, an 
improvement in surface finish should be considered. However, the 
magnitude of the rms roughness and the processes which might best 
improve the surface can not be adequately chosen until the Hertz 
contact diameter is known so a wave length range for measuring surface 
roughness can be defined. 

5 . Wear Measurements by Ferrographic Analysis 

In order to confirm the dependence of asperity interactions 
on A, a ferrographic analysis of the oil used in the test rig was 
performed. In ferrographic analysis, solid particles are precipitated 
from the oil by a magnetic field as the oil flows slowly over a 
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substrate. Particles are precipitated from the oil at a rate 
depending on size and magnetic permeability. The resulting ferrogram 
yields the size distribution and concentrations of particles in the 
oil sample. 

Using special precautions to avoid oil contamination by particles 
other than contact wear debris, 30 m£ samples were taken at various 
Hertz pressures for 1 m/s sliding. The fluid was the same naphthenic 
base oil used in the temperature measurement experiments. Knowing the 
film thickness and surface roughnesses for each of these conditions, 
the ferrogram results could be interpreted in terms of the A parameter. 
Figure III-9 shows the percent area covered on the ferrogram (proportional 
to concentration of particles) as a function of A. Results are 
shown for experiments using the rough, medium rough and smooth balls. 

The two separate curves shown represent positions on the Ferrogram. The 
upper curve corresponds to relatively large size particles compared to 
the second curve. 

The importance of the paramter A is clearly shown. Wear rates 
for both particle size classes (upper and lower curves) begin to increase 
rapidly as A approaches and then falls below one. Therefore, the onset 
of asperity interactions predicted by the occurrence of surface tempera- 
ture fluctuations correlates very well with the results of the Ferro- 
graphic analysis. The calculation of the A value was based on a 
roughness cut-off of 0.76 mm, which is reasonably close to 2d, (0.86 mm), 
the recommended value. 

E. Experiments in Rolling EHD Contacts 

Prior to the investigations being reported here, all EHD 
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experiments conducted in this laboratory were for case of pure 
sliding. The pure sliding experiments were selected because of 
their relative simplicity and because the contact temperatures would 
be higher, thus aiding in the development of the infrared temperature 
measurement technique. However, since a majority of Hertzian contact 
conditions occur in devices having relatively small amounts of slip, 
this condition had to be investigated. 

The apparatus shown in Figure III-2 was constructed to allow 
independent selection of the surface velocity of both the sapphire 
flat and the steel ball. The results presented here are primarily 
concerned with the influence of slip on lubricant film thickness, 
traction and ball surface temperature. 

1. Film Thickness and Traction Measurements 

The results of the film thickness and traction experiments 
may be obtained simultaneously, but they cannot be obtained concurrent 
with the temperature measurements. In addition to requiring two 
different microscopes (metallurgical and infrared) , the film thickness 
measurement requires a partially reflective coating, which would 
interfere with the infrared microscope operation. 

Table III-l is a summary of the film thickness and traction 
experiments. Since only the influence of slip is being studied, a single 
fluid, the naphthenic oil Nl, was investigated. Rolling velocity, 
slide-roll ratio and peak Hertz pressure were all varied. A slide- 
roll ratio I = 0 implies pure rolling and a ratio £ = 2 implies pure 
sliding. In all cases, fluid Nl the naphthenic base oil and the smooth 

ball (.011 um Ra) were used. 
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Film thickness and traction have been plotted as a function of 

E in Figures III-10 and III-ll respectively. The effect of the slide/roll 
ratio on film thickness is small. In fact, both the centerline and 
minimum film thickness depend on sliding velocity to a power in the 
range 0 to -0.09 depending on the Hertz pressure and rolling velocity. 

It is apparent from Figure III-ll that the peak normally obtained in TC 
vs. E experiments in line contacts is not present here. Wedeven's traction 
data for point contacts [23] also lacks such a peak. 

The Influence of rolling velocity on film thickness and traction 
under conditions of zero slip is shown in Figure III-12. 

2. Surface Temperature Measurements 

Table III-2 shows the results of the ball surface temperature 

measurements. Because of small changes in the lubricant bath temperature 

during the course of the experiments, the temperature rise AT (AT = T -T ) 

b bath 

has been tabulated instead of the absolute temperature of the surface. 

Also, due to the relatively low infrared detector signal level under 

conditions of small amounts of sliding, only the high load conditions 
of those listed in Table III-l were imposed in the temperature measurement 

experiments. However, these conditions are the closest to those found 
in normal practice. 

The data from Table III-2 has also been plotted and is shown in 
Figures III— 13— 15. Figures III— 13 and III— 14 show the dependence of 
temperature rise on slide-roll ratio E for rolling velocities of 0.75 
.and 1.00 m/s respectively. A linear regression analysis of the maximum 
temperature rise data for both rolling velocities has resulted in 
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AT = 30.3Z + 1.2, C 


(13) 


with a correlation coefficient of 0.99. 

Figure III-15 shows the temperature rise as a function of 
rolling velocity for the case of pure rolling (£ - 0) . The data fit 
the straight line 

AT = 5 .2V, C (14) 

with a correlation coefficient of 0.99. 

The variation of surface temperature rise with position in the 
EHD contact is shown in Figures III-16-18 for slide roll ratios in the 
range 0 to 2.0. The data for V = .75 m/s (Figure III-16) and V = 1.00 
m/s (Figure III-17) are quite similar and the pure sliding results 
(2 = 2.0) are in good agreement with the data reported previously [13], 
It is interesting to note the small amount of inlet heating present 
under these conditions. For example, at 1.5 Hertz radii upstream of 
the contact center, the surface temperature is less than 5 °C above the 
bath temperature. Considering that the ball surface temperature 
several Hertz radii upstream must be greater than or equal to the bath 
temperature, the 5 °C value is probably only partially due to heating 
in the inlet region. This relatively low value is consistent with the 
film temperature rise in the inlet predicted by Greenwood and Kauzlarich 
[24]. 

A second observation from Figures III-16-17 is the movement of 
the peak in the temperature profile toward the exit as Z is increased. 
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Also, Figure III-18, which contains only zero slip data, shows that the 
peak occurs near the contact center and the temperature at the inlet and 
exit contact boundaries are nearly equal. The trend toward a symmetric 
temperature distribution as I approaches zero is a result of the primary 
mechanism of heat generation going from viscous shear when appreciable 
slip is present, to compression heating at I I 0. Viscous shear results 
in a heat flux present throughout the contact region. If this flux were 
uniform, the ball surface temperature would continue to rise as it passed 
through the contact. The fact that the heat flux is very much greater 
at the contact center results in a temperature reduction as the exit 
boundary is approached. The faster the heat flux falls off, the lower 
the exit temperature should be. In the case of nearly pure rolling, 
the viscous shear component is absent in Hertzian region and the 
temperature rise is due to shear in the inlet and compressional 
heating. The latter term is a function of the pressure profile, which 
is nearly symmetric in the contact. The film and surface temperature 
should therefore rise and fall with pressure. Therefore, the trends 
shown in Figures III-16-18 appear to be reasonable. 

3. Extension of Flash Temperature Analysis for Sliding and 
Rolling 

The flash temperature analysis done in a previous section can 
be extended to consider the influence of rolling. The formulation by 
Archard is followed except that the influence of a second surface must 
be taken into account. Also, since most of the data in Table III-2 
corresponds to L > 5, the fast moving heat source limit is applied. 

It is also assumed that the viscous dissipation takes place in a single 
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plane. This assumption is reasonable because the dissipation is inversely 
proportional to viscosity and the shear stress is constant through the 
film [18]. Since a temperature profile through the film was not 
calculated, it was further assumed that this adiabatic plane occurred 
at the center of the film. If the two surface temperatures are not 
substantially different, this approximation should be reasonable. 

The total heat flux q is the sum of that transferred to each 
surface (assuming convection can be neglected) . Therefore 


TC-W-V 


q - q, + q 

b sa 


s 


TTa 


(15) 


Since the dissipation is assumed to occur only on the midplane, the remain- 
ing fluid merely conducts the heat to the adjoining surface. This allows 
a linear temperature profile to be imposed between the midplane and 
boundary. If T^^ is the midplane film temperature (maximum through 
the thickness) and T^ and T ga are the surface temperatures at the same 
position in the contact, then the conduction equation reduces to 

T - T = 5 (16) 

max sa 2k^ sa 


T - T 
max b 



(17) 


The average surface temperatures can be obtained from the flash temperature 
theory of Archard [18]. It should be noted that the location of viscous 
dissipation is fixed with respect to the IR detector. Therefore, the 
heat source is moving with respect to points on the two surfaces 
with velocities equal in magnitude to the absolute surface velocities. 
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For L > 5 


K 

T - T = q 

sa osa y 1/2 n sa 

sa 


( 18 ) 
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v l/2 q b 
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(19) 


where is the bulk temperature assigned to each surface and 
1/2 1/2 

K = .968a c. /k. Equations 15-19 can be combined to yield AT, = T, - T , . 

b b ob 

The procedure is to eliminate T from equations 16-17 and then to 

max 

substitute expressions for q^ and q^ from equations 18 and 19 into 
the new equation and equation 15. This results in two equations In 
the unknowns, T and T . Eliminating T from them results in the 

S3. D S3 

following equation 



For the conditions used in this study, (T - T . ) is small compared 

osa ob r 

to the first term in the numerator. Using an average film thickness 
of 0.16 ym (Figure III-10) and W = 67N, equation 20 can be evaluated 
for the materials used. 

The ball surface temperature rise given in equation 20 assumes 

only viscous shear as a heat generation mechanism. However, Table III-2 

shows that even when V =0, AT, ^ 0. In order to account for this 

s b 

additional mechanism, the empirical relation for the pure rolling data 
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obtained earlier (equation 14) is used. Adding this temperature rise 
to the rise due to viscous shear (equation 20) results in 


AT, = 5.2V + 


66.7 TC*V 

“^172 ' “ v l/2 

sa _b 

5.8V X ^ 2 +14.9 11,8 

sa 


( 21 ) 


Equation 21 has been evaluated using the conditions shown in Table HI-2 
except for the conditions of pure sliding (E = 0, V =0). The simple 

S cl 

sliding data was excluded because L =0, thus not meeting one of the 

S3. 

imposed restrictions. Also, with a stationary sapphire surface, T^ and 
T ga are likely to be considerably different, resulting in the assumption 
of the adiabatic plane being at the center of the film being in error. 

The results are shown in Figure III-19. 

It should be remembered that the predicted temperature rise is 
the average over the contact area whereas the measured temperature 
rise is the maximum from Figure III-19, it is therefore clear that the 
predicted values of temperature rise are higher than the measured values. 

There are a number of differences in the development of flash 
temperature theory and the conditions found in the experiments. 

First of all, the theory is based on a uniform heat flux over the 
contact area. The heat flux distribution in the EHD contact is probably 
non-uniform with the maximum occurring near the center. Also because 
of differences in thermal dif f usivities of the materials of the contact, 
the adiabatic plane of the fluid film may not occur at the film center. 
Finally, the theory assumes a semiinfinite solid with the free surface 
being adiabatic except under the concentrated heat flux. This is not 
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consistent with the conditions of the experiment. The ball and 
sapphire surfaces outside the contact area are constantly flooded 
with lubricant. The agreement shown in Figure III-19 should there- 
fore be viewed cautiously. 

Because of the apparent agreement shown in Figure III-19, 
equation 20 can be considered as a predictive equation for the 
maximum ball surface temperature. 
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IV. GLASS TRANSITIONS IN LUBRICANTS 
A . Introduction 

A preliminary investigation into the possible role of glass 
transition and glassy state behavior of lubricants in EHD contacts 
was undertaken. Part of this effort requires knowledge of the pressure 
and temperature distributions in the EHD contact which are the subjects 
of the previous sections of this report. In this section preliminary 
results of the glass transition of lubricants as a function of pressure 
by two methods are presented along with a discussion indicating possible 
implications of the results to EHD lubrication. 

The proposition that the mechanical behavior of lubricants in EHD 
contacts might not be that of a viscous liquid but of an elastic solid was 
first presented by Smith [25] in 1960. The results of Johnson and Roberts 
[26], and Johnson and Cameron [27] suggest that the lubricant film responds 
as an elastic solid under the conditions of pressure, temperature and 
shear rate of their EHD experiments. 

Lubricating oils contain linear and highly branched hydrocarbons. 
Thus, a lubricating oil can display both a glass transition which is 
associated with the non-crys tallized (non-linear hydrocarbons) part of 
the lubricant oil and a phase transition which is associated with the wax 
in the oil (linear hydrocarbons) . Regardless of the wax content in a 
given oil, the glass transition temperature remains unchanged [28]. 

When comparing the conditions to which a lubricant is subjected in a 
typical EHD contact to available glass transition data, it appears quite 
likely that the lubricant is in or near the glassy state. This appears 
to be certain for some bulk polymers and polymer solutions and may also 


37 



be true for low molecule weight hydrocarbons. If in or near the glassy 
state, the solidlike behavior and ultimate mechanical properties of the 
lubricant may be the controlling material characteristics in the large 
strain and large strain rate conditions that exist in EHD. 

A preliminary review of the literature [29,30,31] of glassy 
solids indicates that the stress and strains occurring in EHD contacts 
exceed the ultimate stress and strain in many materials. The work per- 
formed in this laboratory on the molecular degradation of lubricants in 
EHD contacts clearly demonstrates that yielding on a molecular level does 
occur in these applications. The infinitesimal strain work of Lamb [32] 
and Litovitz and Dill [33] which is being applied to EHD through the efforts 
of Harrison and Trachman [34] are important in determining the extent to 
which the lubricant is in the glassy state, but, it seems to us, the 
important material properties are the ultimate or large strain properties 
and not the linear viscoelasticity properties. 

In the early stages of this research, glass transition temperatures 
at atmospheric pressure were measured by the technique of Differential 
Scanning Calorimetry (DSC) for two lubricating oils: naphthenic base oil 

R-620-15, designated as N1 and five ring polyphenyl ether, 5P4E. The 
instrument employed was a Perkin-Elmer Model DSC-2. Glass transition 
temperatures were -59 C for N-l sample and -21 C for 5P4E. Probably 
because these values are low compared to the operating temperatures in a 
typical EHD contact, the occurrence of the glass transition phenomenon 
and the existence of the glassy state in the contact were not investigated 
earlier by other researchers in this field. However, because of the fact 
that' the glass transition temperature increases with pressure, and because 
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of the existence of high pressures in the EHD contacts, the glass 
transition temperature can be greater than the operating temperature 
in the contact. The ASME Pressure-Viscosity-Report [35] was used to 
obtain an estimate of the variation of with pressure. The solidifi- 
cation points indicated in the ASME report were assumed to be approximate 
glass transition points. Since these are reported as pressures at 
several different temperatures, they permit an estimate of the change of 
glass transition temperature with pressure. Also, according to McKinny 
and Goldstein [36], the glass transition is an isoviscous state. This is 
apparently because the glass transition is a constant relaxation time 
state for a given process and the relaxation time is essentially 
proportional to the viscosity. The change in temperature with pressure 
required to maintain constant viscosity gives results that are consistent 
with the results based on solidification. Based on this approach for 
typical lubricants, the glass transition temperature increases with 
pressure at a rate ranging from 80 to 350 C/GPa. Thus, it would be 
expected that some lubricants will be in the glassy state in an EHD contact 
with average pressures of 0.7 GPa or higher at room temperature. There- 
fore, it would be expected that many lubricants are in the glassy state 
for a significant portion of the time they are in the contact. 

The temperatures and pressures at which glass transition occurs 
and the mechanical properties of the lubricant in the glassy state may 
influence the two most important dependent operating variables of an 
EHD contact: the film thickness and the traction. Since the relaxation 
time of a material in the glassy state is long compared to the residence 
time in the contact, it is possible that once in the glassy state, the 
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lubricant will remain in the glassy state on the moving surface while 
the surface moves from one EHD contact to the next. This may happen 
in spite of the fact that the pressure causing the glassy state forma- 
tion has been removed. It is also possible that the glassy lubricant 
does not have time to flow while in the EHD contact but cleaves some- 
where in the film so that solid layers of lubricant adhere to the contact 
surfaces and slide upon each other or the lubricant may undergo large 
strain like a ductile solid subjected to large stresses. 

Therefore, the glass transition temperature as a function of 
pressure is a significant material property of lubricants. The 
mechanical and thermal properties of a substance in the glassy state 

are known to be functions of T (P) which in turn is a function of 

g 

temperature, pressure, and rate of temperature and pressure changes 
prior to entering the glassy state. 

B. Glass Transition Temperature and the Glassy State 

Glass transition is characterized by certain experimental 
observations which occur while the imposed environment of the material 
is changing. Commonly experiments are performed for isobar ic 
cooling or isothermal compression of the material. At the glass transi- 
tion, many material properties change in a characteristic manner. As the 
temperature is reduced in the liquid region, the material contracts and 
the viscosity increases. If the material is capable of crystallization, 
a point will be reached where crystallization starts. If crystallization 

does not occur or can be avoided, the viscosity will continue to increase 

7 12 

until some level of 10 -10 Pas is reached [37,38,39]. At this point, 
the material becomes rigid and the thermal expansion coefficient falls by 
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about one-half to one-third of its value in the melt or liquid state 
as shown in Figure IV— 1. 

Glass transition is also accompanied by a change in the specific 
heat, 0 ^, isothermal compressibility, 3, and other secondary properties. 

It corresponds to a change in slope of a plot of specific volume versus 
temperature as shown in Figure IV-2. Thus, glass transition phenomena 
are a characteristic of any liquid which can be supercooled to a suffi- 
ciently low temperature without crystallization. 

Glass transition phenomena are often referred to as an apparent 
second-order transition since it is characterized by a discontinuous 
change in the secondary thermodynamic quantities. These changes occur 

over a range of temperatures and are not strictly discontinuous. 

Ferry [403 and Haward [29] showed that glass transition phenomena are 
not true second— order thermodynamic transitions since at T^ the substance 
is not in thermodynamic equilibrium. This is due to the slowness of 
molecular rearrangements at this temperature. This absence of thermo- 
dynamic equilibrium is part of the definition of the glassy state. 

Below T , the degree of order will appear fixed and will not 
vary with temperature or pressure during the time of experimental 
observation. Because of this, the structural degrees of freedom are 
said to be frozen-in and therefore, the structural contribution to a, 

c and 3 are absent in the glassy state. 

P 

The effect of pressure on T is shown in Figure IV-2. The glass 

S 

transition temperature shifts upward as the pressure is increased. The 
increase in T^ with pressure will be sufficient to maintain an isoviscous 
state. A sufficiently high pressure can induce a transition from the 
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liquid state to the glassy state without requiring a decrease in the 
temperature. 

The transition temperature of a material is a function of the 
imposed rate of change of the material's environment, and for a series 
of constant rate experiments, the transition temperature curve is a 
constant relaxation time curve. The characteristic time of the experi- 
ment is dependent on the observation time employed in the experiment. 

An approximate value of the relaxation time of the material when 
transition occurs is the observation time employed in the experiment. 

The range in the relaxation time of the material as it goes through 

2 

the transition from liquid to glass is from about 0.1 to 10 times the 
observation time. In the transition region the material exhibits both 
viscous and elastic behavior relative to the process. The range of the 
transition region is indicated on Figure IV-1 and 2. 

As a result of the non-equilibrium state, the thermodynamic 
history of a glass forming liquid has considerable influence on the 
structure, the transition and the properties in the glassy state. The 
rate of the imposed environmental change effects the transition point and 
the properties of the glassy substance. Glass transition can be reached 
by isobar ic cooling, isothermal compression, an imposed rate change or 
a combination of these changes on the material. For example, the 
influence of the rate of cooling, K, on the shift of the glass transition 
temperature of poly (vinyl acetate) has been studied over a wide range of 
rates by Kovacs [41]. The volume-temperature response obtained from 
different constant rates of cooling of a glass-forming substance is shown 
in Figure IV-3. Since decreasing the rate of cooling increases the 
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effective experimental time, the glass transition temperature will be 
shifted to lower values. Matsouka and Maxwell [42] studied the effect 
of rate of compression on the compressibility curves of polysterene at 
121 C. Their results indicate that glass transition pressure, P , 
shifts to lower values when a greater rate of pressure application is 
employed. If the pressure is increased slowly, the molecules of the 
material will have enough time to rearrange and the glass transition 
will take place at a higher pressure. On the other hand, if the rate 
of pressurization is increased, the time for any structural changes is 
smaller and the glass transition occurs at a lower pressure. 

The influence of thermodynamic history on glass transition of 
poly (vinyl acetate) was studied by McKinney and Goldstein [36] by using 
three different thermodynamic histories: variable formation history 

(isobaric cooling at different pressures) , and two constant formation 
histories (cooling at atmospheric pressure or 80 MPa (800 bar) followed 
by pressure changes in the glassy state) . 

Several techniques are employed for the determination of the 
glass transition temperature such as dilatoraetry [36,38,43], differential 
scanning calorimetry [28,44,45], thermomechanical analysis (TMA) , 
dielectric [46,47] and light-scattering techniques [48-58]. Depending 
on the technique used, various methods exist for specifying the glass 
transition temperature within the transition region. The usual dila- 
tometric technique is to cool the liquid at a constant rate and 
extrapolate the linear portions of the volume-temperature relation 
above and below the transition region to their intersection. The 
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temperature of this intersection [36] is taken as the glass transition 

temperature as shown in Figure IV-2. If T is measured by the change in 

g 

the expansion coefficient or the specific heat, then is taken as 
the mid-point in the step-change as measured from the extensions of the 
glass and liquid base lines as shown in Figure IV-1. While the choice 
is somewhat arbitrary, and other authors have suggested alternate tech- 
niques, the above methods are the most commonly used. Because of these 
different ways of defining T for different experiments, and because of 

o 

the different inherent rates and histories in the different experiments, 
different glass transition temperatures may be measured on the same 
material in different experiments. 

The first step in the study of the importance of the glassy 
state in EHD contacts is to determine if the lubricant is in or near 
the glassy state while in an EHD contact. This requires the determina- 
tion of the glass transition temperature as a function of pressure and 
rate at observation times comparable to that of an EHD contact for typical 
lubricants. From these data, a phase diagram showing the liquid and the 
glassy states for the lubricant can be constructed. Based on the tempera- 
tures and pressures in the contact inferences can then be made regarding 
the state of the lubricant in the contact. 

Knowing the temperature and pressure distributions in the EHD 
contact, the phase diagram will permit determination of whether or not 
the lubricant is in the glassy state, approximately what operating 
conditions result in the glassy state behavior, and over what fraction- 
of the contact area the glassy state exists. 
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C . E xperimental Technique 


The glass transition temperature as a function of pressure based 
on two different techniques - light scattering and volume dilatoxnetry - 
is reported. The two techniques are complimentary because the observa- 
tion times are similar (about 200s) but the characteristic rates of the 
experiments are quite different - (10 + ^s ^ for dilatometry and 10 ^s ■*- 
for light scattering) . Therefore as the material approaches the amorphous 
glassy state in the case of the volumetric experiment the material is a 
viscous liquid while in the light scattering it is behaving as an elastic 
liquid . 

1. Light-Scattering Method 

The basic concepts of the light-scattering technique and the 
description of experimental equipment used for measuring the glass 
transition temperature are discussed. 

a. Theory 

When a laser beam passes through a transparent medium, a small 
portion of the light is scattered in all directions. This light 
scattering is a result of optical inhomogeneities in the scattering 
medium. The optical inhomogeneities are the result of fluctuations of 
the dielectric constant brought about by fluctuations of density 
due to random thermal motion of the molecules [59] . 

The dielectric constant, e, of a liquid can be considered as 

e = e + Ae (1) 

where e is the average dielectric constant and Ae represents the 
fluctuation in e. Assuming the entropy, S, and the pressure, P, are 
independent variables, one can write Ae as a sum of two terms 
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( 2 ) 


ie = (ff) p AS + f|f) g iP 

The first term is the result of entropy fluctuations at constant pressure 
(isobaric) while the second term is the result of pressure fluctuations 
at constant entropy (adiabatic) . The convenience of this choice of 
independent variables is obvious if one realizes that thermal sound 
waves are adiabatic pressure fluctuations. These pressure fluctuations 
represent random local compressions or rarefections which, as a 
consequence of the elastic properties of the medium, do not remain fixed 
in position but travel throughout the volume of the sample. These 
adiabatic fluctuations can be described by means of plane sound waves 
of thermal origin propagating in all possible directions. 

The study of the specific heat of a solid lead Debye [60] to the 
conclusion that the energy which one must have in all 3N degrees of 
freedom associated with the atomic oscillations of the solid can be 
considered as the energy of 3N normal elastic waves. Thus, Debye 
considered the energy of thermal motion of a solid as the energy of 
elastic waves. The Debye representation of the energy of thermal motion 
can also be extended to liquids and dense gases if the damping of thermal 
waves is small. 

A major advancement in the development of light-scattering theory 
was achieved when Brillouin [61] established the fact that sound waves 
and the Debye thermal elastic waves were the same phenomena. 

Thermal fluctuations do not remain "frozen", but change continu- 
ously in time and therefore, they can modulate the scattered light. If 
we let the modulation function be $(t) and the monochromatic light wave 
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undergoing modulation be given by 

E (t) = E exp(iu) t) (3) 

s o o 

where E (t) is the electric field intensity of the scattered light and 
s 

E is its amplitude, the resulting vibration of the scattered light 
o 

will have the form 


E’(t) = E (t) • $>(t) (4) 

s s 

Consider the modulation of the scattered light by adiabatic fluctuations 
of pressure. The change of pressure with time satisfies the wave 
equation. Neglecting the absorption of the elastic waves, we can write 
the wave equation in the form [59] 

$ - V 2 V 2 $ = 0 (5) 


where V represent the velocity of propagation of the elastic waves. 
The solution of equation (5) is given by 


0(t) = cos [ftt - k- r ] 


( 6 ) 


where ft is the frequency of the elastic wave and equal to V*k where k 
is the wave vector of the elastic wave. Thus, the vibration in the 

-f 

scattered light wave with wave vector k g diffracted by the elastic wave 

4 - ■ 4 ' 

with wave vector k in the direction of k g is given by 


E 

E^(t) = — {exp[i(u) o + ft)t - k*r] + exp[i(a} Q - ft)t + k*r]}(7) 


The above result shows that there are two doublets in the spectrum 

of the scattered light with frequencies (oj + ft) and (u) - ft) 

o o 
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symmetrically located on the two siaes ui uie incident frequency, 

All 3N Debye waves in a body propagate in all possible directions. 
If such a body is irradiated by a parallel light beam with wave vector 
and the scattered light is observed in the direction of the wave vector 
of the scattered light, k g , making an angle 0 with k^ as shown in Figure 
IV-4, then the maximum intensity of the scattered light along k g will 
occur when the wave vector k of the elastic wave together with the vectors 

— h 

and k^ satisfy the Bragg condition 

-> 

k. + k = k 
1 s 


i ->■ — —y ->■ 0 

k = k = + k - k. = 2n„k. sin — 
1 s l i l 2 


where n^ is the refractive index of the scattering medium at the incident 
wave length. Since the wave number is the reciprocal of the wave length, 


equation (8) can be rewritten as 


2n , sin 0/2 £1 K J 

i 

where A is the wave length of the sound wave and X is the wave length 
of the incident light, while V is the sound velocity and fi the frequency 
of the sound wave. 

Moreover, as a consequence of the sound wave, the frequency of 
the incident light is altered by the Doppler effect with respect to 
the original value oj . This shift was first predicted by Brillouin [61] as 


= + 2o) f sin 0/2 
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where c Is the speed of light. 

Rearranging and using the relation X = c/to^ we can write equation 

(10) as 


V = + OL. 

B 


A 

2n^sin(0/2) 


(ID 


Comparing equations (9) and (11) we obtain the frequency shift, u> , equal 

to the sound frequency, fi. For scattering angle of 90° and incident light 

of wavelength of 5145°A, the frequency shift in materials having a sound 

velocity of the order of 3000 m/s. and a refractive index of 1.5 is about 
9 

12 x 10 Hz. This result shows that the Brillouin method deals with 

acoustical phenomena in the gigacycle region (hypersound) . This shift 

represents a frequency change of less than 0.002 percent. The 

observation of this frequency shift is possible with a gas laser as a 

monochromatic light source and a Fabry-Perot interferometer as a 

detector. The two signs in the above equations correspond to the two 

opposite directions in which the sound wave may travel. 

According to Landau and Placzek [62], the scattered light 

produced by the entropy fluctuations of Ac does not propagate in normal 

liquids and these fluctuations are the source of the central unshifted 

component of the scattered light. The intensity of this component, I , 

c 

is expressed as a function of the total intensity in the Brillouin 

doublet, 21 , by the Landua-Placzek relation [48,49] 

B 



Y-l 


where y is the specific heat ratio. 


( 12 ) 
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Measurements on Brillouin spectra have been made by a number of 
workers such as Rank, et al. [48,50], Cummins and Gammons [51], Pinnow 
et al. [49], and others [52,53] who observed the Brillouin effect by 
using He-Ne and Argon ion laser as the light source. A typical frequency 
spectrum of lubricating oils is shown in Figure IV-5. The shift of the 
Brillouin lines from the laser frequency is of interest since it enables 
us to calculate the velocity of sound according to equation (11) . 

For a viscous medium,, two types of Brillouin scattering can be 
observed ' 

(i) Scattering from longitudinal sound waves which gives a 
measurement of the longitudinal modulus, M, as 

M = pvj[ (13) 

where is the longitudinal velocity of sound given by 

W BJ?, 2n.sin 0/2 (14) 

l 

(ii) In highly viscous materials just as in crystals and glasses, 

a second set of Brillouin lines is often observed in the depolarized 

spectrum with a shift smaller than that of the longitudinal Brillouin 

lines (not shown in Figure IV-5) . These new shifted components are 

transverse Brillouin lines due to scattering from propagating thermal 

transverse waves which can propagate in the medium as it becomes more 

viscous. Scattering from transverse waves appear in the depolarized 

spectrum because these waves cause anisotropic fluctuations in the local 

structure of the material. In lubricants at high viscosities, the 

material is acting in a glassy manner and the transverse Brillouin 
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peaks may be detected. The fact that these peaks appear in the frequency 
spectrum is an indication that the longitudinal waves are no longer 
purely compress ional waves but contain a transverse component also. 
Scattering from transverse sound waves gives a measurement of the shear 
modulus, G, as 

G = p V* (15) 

where V is the transverse velocity of sound given by 

V t = W Bt 2n.sin 0/2 (16) 

1 

In low-viscosity liquids, the transverse elastic waves are so weak that 
the inner Brillouin doublet, which is brought about by the diffraction 
and modulation of light by the transverse wave, is difficult to observe. 
In this research only longitudinal frequency shifts and longitudinal 
sound velocities will be reported. 

b . Pressure and Temperature Coefficients of the Sound Velocity 
A change in the slope of the sound velocity as a function of 
temperature and pressure like density or specific heat also defines the 

glass transition and has been observed by many workers [52-55]. It has 

been established that the change in the temperature coefficient of the 
sound velocity is a manifestation of the change in the thermal expansion 
coefficient that occurs at the glass transition temperature [56-58] . 

This allows us to detect the glass transition with gigacycle sound waves. 
Figure IV— 6 shows the variation of the sound velocity and frequency shift 
with temperature for atactic polystyrene [53] where a change in slope 
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is clearly observed at the glass transition temperature. 

The dependence of the velocity of sound on the density, p, is 
seen through the relation 

V = (M/p) 1/2 (17) 

where M represents the appropriate modulus of elasticity which is a 

function of the density. Since the density appears explicitly in the 

above relation, it is clear that a change must appear in the slope of 

the velocity-temperature and velocity-pressure curves at the glass 

transition. A comparison of volume-temperature data with velocity- 

temperature data taken from literature [56] is consistent with the above 

observation. However, the expansion coefficient of polystyrene is 

2.2 X 10“ 4 V 1 for T < T and 5 x 10 _4 °C~ 1 for T > T while the 

g g 

—4 o —1 —4 

velocity coefficient is 8 x 10 C for T < T and 14 x 10 for T > T 

g g 

which shows that the dependence of velocity on density is chiefly through 
the strong dependence of the modulus on density and only to a lesser extent 
through the explicit appearance of density in equation (17) . 
c. Light Scattering Experimental Apparatus * 

The schematic arrangement of the basic components of the light- 
scattering experiment is shown in Figure IV-7 and consist of (1) a 
argon ion laser in single frequency operation and having a 
power of about 200 milliwatt at X = 5145°A (single line operation is 
accomplished by inserting a wavelength-dispersing element (Prism) into 
the laser cavity), (2) a high pressure scattering cell containing the 

*The assistance of Drs. D. C. O’Shea and F. D. Medina of the 
Georgia Tech Physics Department for the use of their light scattering 
equipment and guidance in this work is gratefully appreciated. 
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lubricant sample to be investigated, and (3) a detector which permits 
the spectroscopic analysis of the scattered light. 

The scattered light from the lubricant sample is collected by 
a lens and enters a Fabry-Perot interferometer as shown in Figure IV-8. 

The output of the interferometer, consisting of polarized and depolarized 
scattered light, is focused on a screen which has a pinhole in front of 
a photomultiplier tube (PMT) . The output signal from the PMT can be 
processed by a standard counting system and recorded with a multichannel 
analyser. The spectrum can also be recorded on X-Y recorder and/or 
punched on paper tape using the output of a rat erne ter. 

The Fabry-Perot interferometer used in the light-scattering 
experiment consists of a pair of highly reflective mirrors separated 
by a specified distance, d. The mirrors can be moved and tilted to 
ensure exact parallism between their surfaces. At certain wavelength [63] 
the amount of light transmitted through the mirrors will be many times 
larger than expected from the combined transmission for two highly 
reflective mirrors. This occurs when there is constructive interference 
between incident light waves with those of different reflections. Only 
when all the waves from the multiple reflections are in phase will there 
be high transmission. For incident light normal to the mirror faces, 
the maximum transmission will occur when 



m m 


(18) 


where m is the order of interference (an integer). In terms of frequency, 
the above equation can be written as 
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m 


( 19 ) 


The transmitted beams are focused by a lens onto a screen where a series 
of concentric interference fringes are observed. Only light travelling 
perpendicular to the plates obey the simple resonance condition. This 
light is focused to the center of the ring pattern. By placing a pin- 

Q 

hole at this center, only light with frequency V = m — is allowed to 

m Za 

pass through and be detected by a photomultiplier tube. If we plot the 
transmission of the interferometer as a function of frequency for a 
monochromatic incident light of frequency f, we would obtain a curve 
similar to that shown in Figure IV-9. These transmission characteristics 
enable the device to be used as a tunable filter to display the output 
frequency of the laser. The frequency range which is transmitted during 
the Fabry-Perot scan is known as the free spectral range (FSR) which is 


FSR = (20) 

O 

Since the speed of light is approximated by 3 x 10 m/s, the frequency 

units can be expressed in cm ^ (wave number units) . A frequency of 
9 -1 

15 x 10 Hz corresponds to 0.5 cm . Unless otherwise specified, the 
frequency units will be reported in cm ^ . 

In the scanning mode, we observe the spectral content of the 
scattered light by moving one mirror along the interferometer axis with- 
out disturbing the parallelism of the mirrors by applying a ramp 
(sawtooth) voltage to a piezoelectric transducer (PZT) attached to one 
of the mirrors . 

The temperature inside the Fabry-Perot interferometer is controlled 
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by a proportional controller and a thermister sensor. Two resistors 
(20 watt) are glued to each of the bottom and side panels of a box 
enclosing the interferometer and serve as heaters. The temperature is 
controlled to about 2-3° above room temperature. 

Since the Rayleigh line intensity is usually several orders 
of magnitude greater than that of Brillouin lines, a significant improve- 
ment of the frequency spectrums can be obtained by operating the Fabry- 
Perot interferometer in triple pass geometry. This is accomplished by 
using a pair of corner cubes (Figure IV-8) . 

In summary, the light-scattering technique enables one to 
determine the sound velocity at frequency if a monochromatic beam of 

D 

light strikes the scattering medium and the frequency shift of the 
Brillouin scattered light at angle 9 is measured. The location of the 
change in the slope of the sound velocity as a function of temperature 
and pressure is a measure of the glass transition. 

d. High Pressure Equipment 

The high pressure equipment consists of a scattering cell, a 
strain gauge pressure transducer, a 9:1 area ratio intensifier and a 
hand pump . 

The high pressure cell (3" x 2-3/4” x 1-3/4”) is shown schematically 
in Figure IV-10. It was used to pressures of 0.66 GPa (100,000 psi) . The 
test fluid in the cell was separated from the hand pump oil by an intensi- 
fier which was connected to the cell by high pressure tubing. The cell 
material is Vasco Max 300 CVM alloy. 

*The high pressure cell was originally designed by Dr. W. B. 

Daniels of the University of Delaware. 
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Heat treatment of the cell was necessary to raise the hardness 
and increase the strength properties. Aging the cell by heat treatment 
for six hours at 480 C followed by an air quench raised the hardness to 
about 52-55 Rockwell "C" and increased the yield strength to 1.9 GPa. 

An intensifier with a 9:1 piston area ratio was mounted between 
the hand pump and the high pressure scattering cell. The test fluid is 
contained in the high pressure side of the intensifier cavity, the high 
pressure tubing and the optical cell. The high and low pressure sides 
of the intensifier piston were sealed with Viton 0-rings which were pro- 
tected by anti-extrusion rings. The intensifier material is the same as 
that of the optical cell and it was heat treated in the same manner. 

The pressure of the lubricant sample contained in the optical cell 
was measured with a commercial strain gage pressure transducer located as 
shown in Figure IV-10. The pressure transducer is a Norwood model 114 
manufactured by the Advanced Technology Division of American Standard and 
was used in high pressure viscosity measurements conducted under this grant. 

One of the important features of using the pressure transducer is 
that the lubricant sample pressure is measured directly in the optical cell 
and thus the influence of intensifier seal friction on pressure measurements 
is eliminated. However, when the material viscosity increases during glass 
formation the sample pressure may be different from the transducer pressure 
because of glass formation in the tube connecting them. 

e. Glass Formation History 

Because of the dependence of the glass properties on formation 
history, two standard procedures were adopted to form the glass of the 
materials investigated. In one formation history (history A), the 


pressure was Increased from atmospheric to a reference pressure in the 
liquid region. At the same time, the sample temperature was increased 
at a rate of 6.1 mC/s (22C per hour) to about 17 C above its assumed 
glass transition temperature. At this point the pressure and temperature 
were kept constant for about 30 minutes. The sample was then cooled 
at constant pressure and frequency spectra were taken at 2.8C (5 °F) 
intervals. Twenty minutes were allowed for the temperature to reach 
equilibrium before each spectrum was recorded. Hence, experimental 
measurements were started in the liquid state and the glassy state 
was reached by decreasing the temperature at constant pressure as 
shown in Figure IV- 11. This procedure was repeated at different constant 
pressures and, therefore, the structure of the glass formed varies with 
the formation pressure. 

In the second formation history (history B) , the glass was 
formed by pressurizing the sample at room temperature from atmospheric 
pressure until the glassy state was reached. Frequency spectra for 
this case were taken at pressure intervals of about 20 to 27 MPa 
(3000 to 4000 psi). Again up to twenty minutes was allowed for 
the sample to reach equilibrium before the spectra were recorded. 

In history A, the sample temperature was changed by controlling 
the voltage input to an electric tape heater wrapped around the 
intensif ier— cell assembly which was enclosed in an oven. By controlling 
the input voltage to the electric heater, the cell temperature was 
controlled to better than + 0.1 C. The temperature was measured with a 
copper-constantan thermocouple in conjunction with a direct reading 
digital thermometer. The thermocouple was in a brass well fitted into 
a hole in the cell body (Figure IV-10) . 
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f. Data Reduction Technique and The Re sults for 5P4E 
The data on 5P4E will be used to illustrate the data reduction 
technique. Some results are also available on Nl* and the phase 
diagram for that material based on light-scattering will be presented 
in the Discussion section. 

All frequency spectra of scattered light were obtained at a fixed 
scattering angle of 90°. The free spectral range for each case investi- 
gated was kept constant. This procedure was adopted to minimize the 
errors in measuring the separation distance of the Fabry-Perot inter- 
ferometer which in turn minimized the error in the calculated frequency 
shift. Using the scattering angle of 0 = 90° in equation (9) and (11), 
one obtains 

A = — (21) 

/l n. 

l 

and 


V 


= UL 


X 

/in. 

i 


( 22 ) 


The recorded spectra were used to measure the frequency shift and to 
calculate the sound velocity by employing equation (22) . 

The frequency spectrum of 5P4E recorded at 0.2 GPa (28,200 psi) 
and 24.4 C (76 F) (formation history B) is shown in Figure IV-12. It 
contains two successive orders of the interference. The center of the 
four Brillouin peaks were located at some height where the Brillouin 

*See Appendix A for material descriptions. 
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components appeared to be symmetrical. The separation distance, d, of 
the Fabry-Perot interferometer was 0.50 cm corresponding to a free 
spectral range of 30 GHZ . The frequency shift obtained from this 

spectrum is equal to 


£ . Jr_ 

W B,k 2L, 2d 
k 


(23) 


where k = 1,2. Two values of uj can be calculated, each corresponding 

B 

to the order of interference. These are (see Figure IV-12) 


- & . _c_ 

^B , 1 “ 2L 1 ’ 2d 


6.85 
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3 x 10 
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= 12.45 GHZ 


and 
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6.85 


3 x 10 


10 


16.40 A 2 x 0.50 


- 12.53 GHZ 


Each value of co together with the values of n, (1.6873) and A(5145 A) 
B i 

is used in equation (22) to determine the velocity of sound at the 
indicated pressure and temperature as follow 


Aco. 


v i ■ 


Li = 

/l n . 

1 


-in t C in 9 5145 x 10 
= 12.45 x 10 x 


-10 


1.6873 /2 


= 2685 m/s 


and 


Aco 


V 2 = 


/2 n. 


-men 1^9 5145 X 10 

= 12.53 x 10 x 


-10 


1.6873 / 2 


= 2702 m/s 
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The average values of u) 1 and V are reported. Frequency shifts and 

B 9 k 

sound velocities at other pressures and temperatures are calculated in 
the same manner. The results of these calculations for 5P4E at 24.4 C 
(76 F) (history B) are shown in Table IV-1 and plotted in Figure IV-13. 
From this figure, it is seen that the sound velocity and the frequency 
shift increase with increasing pressure and a change in the slope is 
apparent at 0.17 GPa (24,700 psi) corresponding to a velocity of sound 
of 2658 m/s* This change of slope represents the glass transition. The 
velocities on each side of the transition region can be represented by 
linear functions of pressure. The two sections of Figure IV-13 were 
each fit by least squares regression to a linear expression of the form 

V = A + B P (24) 

where V is in m/s and P in MPa. The intersection of these lines are 
interpreted [52-58] as occurring at the glass transition pressure. The 
least square expression obtained for each section of the velocity- 
pressure curve shown in Figure IV-13 is 

V T = 2098 + 156.3 P 

JLi 

and 

V = 2397 + 72.87 P 
g 

Solving these two equations for their intersection, one obtains 

P = 170.3 MPa and V = 2658 m/s. 
g 

Frequency spectrum of 5P4E obtained by formation history A 
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(at 0.40 GPa (59,000 psi)) and recorded at 64.4 C (148 F) is reproduced 
in Figure IV-14. It is typical in general appearance of all frequency 
spectra obtained. The separation distance, d, was 0.725 cm corresponding 
to FSR of 20.7 GHZ. The frequency shifts and the sound velocities were 
calculated from the position of the Brillouin peaks as discussed above. 
The result of these calculations are shown in Table IV-2 and are plotted 
in Figure IV-15. The sound velocity increases as the temperature is 
decreased. At the transition region the dependence of sound velocity on 
temperature changes. This change of slope of the temperature 
coefficient of the sound velocity was observed at all formation 
pressures. Once again, the two portions of the sound velocity curves 
were fit by least square regression to straight lines of the form 

V = a + b T (25) 

where V in m/s and T in C. 

5P4E was subjected to formation history A at a total of six 
formation pressures. The resulting material behavior is shown in Table 
IV-3 and Figures IV-16 and 17. 

The least square expressions for both portions of the velocity- 
temperature curves were solved for their intersection to determine the 
glass transition temperature at each formation pressure. The glass 
formation temperature-pressure combinations result in the phase diagram 
shown in Figure IV-18. The glass transition temperatures and pressures 
were least square fit with a straight line. The expression obtained 
is given by 
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l 


T = 0.183 P - 4.74 
g g 


(T = 2.27 P +23.5 with T in F and P in kpsi) 
g g g g 

with T in C and P in MPa. The glass transition temperature at 
g g 

24.4 C (76 F) by history B was excluded from the fit due to the different 
history by which the glass was formed. However, in general, the transi- 
tion obtained by History B fell near the line obtained by History A. 

The scattered light frequency shift was also least square 
fit with straight lines and the expressions obtained were solved 
for their intersection to locate the glass transition. The 
results obtained from the frequency shift did not differ from those 
obtained from the sound velocity by more than 1.4 MPa (200 psi) for the 
glass transition pressure (History B, at 24.4 C (76 F)) nbr more than 
1 C (2 F) transition temperature at any formation pressure for History A. 
Therefore, the glass transition data were obtained from the sound 
velocity only. 

D. Dilatometry Method 

As mentioned in the introduction, dilatometry can.be used to 
measure the glass transition at long observation times. In many respects 
volume dilatometry is both conceptually and experimentally less compli- 
cated than the light scattering method for determining glass transition. 
Dilatometry is inherently a low rate process and a change in rate 
shifts the glass transition in a known manner. Therefore, dilatometry 
data is significant in that any temperature-pressure combination 
placing the material in the glassy-state by dilatometry measure- 
ments will indicate glass formation in an EHD contact because of the 
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shorter observation time in the EHD contact. 


Dilatometry measures the relative volume of the material as a 
function of temperature and pressure. The compressibility and thermal 
expansion coefficients of the material differ in the liquid and glass 
regions and the glass transition is defined by this change as discussed 
in the introduction and shown in Figures IV-1, 2, and 3. 

Three dilatometers were constructed and are shown schematically 

in Figure IV-19, IV- 20, and IV-21. The three devices were made to cover 

different pressure ranges. Two types of experiments can be conducted: 

1) isobaric cooling (History A) and 2) isothermal compression (History B) . 

Both the isobaric cooling and isothermal compression experiments had 

2 3 

observation times on the order of 10 s to 10 s. 

The formation histories in the dilatometry experiments were similar 
in that the material was always taken from the liquid state into the 
glassy state. For elevated pressure experiments this required heating 
and pressurizing in the liquid region well above the glass transition to 
some starting pressure and temperature from which the isobaric cooling 
or isothermal compression would start. As the experiment progressed the 
relative specific volume would be obtained as a function of temperature 
or pressure and the intersection of the curves in the liquid and glassy 
regions would define the glass transition. These intersections and resulting 
phase diagrams were obtained in the same manner as that in the section on 
light-scattering for the sound velocity data. 

Experiments have been completed on the materials shown in 
Table IV-4. The phase diagram for 5P4E obtained with the volumetric 
method compared with that from the light-scattering is shown in 
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Figure IV-22. The reason for the difference in the data from the two 
techniques is not known at this time but the difference is small 
relative to the temperature and pressure ranges a material is subjected 
to in an EHD contact. Therefore, if 5P4E is subjected to a temperature- 
pressure combination in an EHD contact which is to the right of the 
dilatometry curve on this phase diagram, the material will be in the 
glassy state in the EHD contact. These data suggest that 5P4E would 
be in the glassy state in most EHD applications. Figure IV-23 contains 
similar comparisons for material Nl. The same comments made regarding 
5P4E are applicable to Nl. 

Figure IV-24 contains phase diagram data for several materials 

based on the isothermal compression measurements. The variation in 

glass forming tendencies of the materials is large. However, the 

typical mineral oil based lubricants studied have very similar glass 

forming characteristics as shown in Figure IV-25. The transition 

regions for the compression experiments are shown in Figures IV-26, 

IV-27, and IV-28. Because the observation time in these experiments 

was approximately five minutes the glass transition point is where the 

material relaxation time is approximately five minutes and the transition 

region represents the changes in relaxation time (x ) from about one 

2 

tenth that value on the liquid end to about 10 times that on the glassy 
end (0.5 <_ x^ j< 500 minutes) . In this region the material exhibits both 
viscous and elastic behavior relative to this experiment. In the 
liquid region the relaxation time is much less than 0.5 minute and the 
material can be treated as a viscous liquid in this process. In the 
glassy region the relaxiation time is greater than 50 minutes and the 
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material can be treated as an elastic solid in this process. 

Table IV-5 lists those materials which were studied but for which 
no glass transition was observed for the volumetric experiments in the 
pressure and temperature range studied. This implies that the relaxation 
times of these materials were much less than one half minute in the 
pressure and temperature range of observation. 

E. Two Additional Experiments 

Glass transition is accompanied by changes in all the material’s 

physical properties including very large changes in viscosity. In the 

past glass transition has been defined as the point at which the viscosity 

11 2 

of the material reaches a high value (e.g. 10 Ns/m ). Because the 
approximate shear rates and shear stresses to which the material is 
subjected in an EHD contact are not consistent with such high viscosities, 
while the temperatures and pressures indicate it is in the glassy state, 
two simple shearing experiments were conducted. Both were conducted in a 
concentric cylinder device. In one experiment the torque (shear stress) 
was measured continuously at constant shear rate (0.43 s while the 
5P4E sample was cooled continuously. In the second experiment a crude 
indication of the ultimate shear strength of the material in the glassy 
state was obtained. 

The viscosity measurement experiment was conducted in a Haake concen- 
tric cylinder viscometer and the results are shown in Figure IV-29. 

As expected the viscosity increases rapidly with decreasing temperature 
but the trend begins to change substantially as the glass transition 
temperature is approached. The glass transition temperature as measured 
by differential scanning calorimetry is -23 C and is marked on the 
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figure. The shear stress limit of the instrument is about 3.5 kPa. As the 
temperature continues to decrease below the glass transition temperature 
the material effectively becomes a rigid solid. 

Once in the glassy state the material has solid-like properties 
such as sound velocity, thermal expansion coefficient, specific heat, 
etc., and might also be expected to have an ultimate yield shear 
strength. To get an estimate of this property at atmospheric pressure, 

a concentric cylinder device was constructed with a outer diameter 
of 16 mm, inner diameter of 13 mm and length of 16 mm. The outer cylinder 
could be held firm and the inner cylinder turned with a torque wrench 
when the sample was cooled below its glass transition temperature. A 
thermocouple was inserted into the sample through the inner cylinder 
and the device cooled by liquid nitrogen. The torque required to break 
the inner cylinder loose was then an estimate of the yield shear strength 
of the material. The data are shown in Table IV-6. These results suggest 
that yield behavior of the material in the glassy state should be studied 
in conjunction with the study of material behavior in EHD contacts. 

F. Discussion of Results 

1. Pressure Dependence of Glass Tr ansition Temperatures in Lubricants 

The glass transition results presented suggest that glass transi- 
tion phenomena and rheological behavior in the glassy state are important 
for some lubricants in some EHD applications. 

The volumetric experiments covered the widest range of pressure 
and temperature. As seen in Figure IV-23 for N1 at low formation pressures 
(up to about 0.45 GPa), the transition temperature is a linear function 
of pressure with a slope of about 190 C/GPa. Above 0.45 GPa, the slope 
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decreased gradually and the glass transition relation became linear 
again above 0.55 GPa with a slope of about 100 C/GPa. Figure IV-23 
also Includes the light scattering transition data obtained by history A 
for fluid N1 (slope = 120 C/GPa) . The glass transition temperatures of 
the light scattering technique are seen to be higher (by about 19 C) 
than the temperatures based on the volumetric method at all pressures. 

The reason for this difference is unclear at present. 

The glass transition temperature for N1 at atmospheric pressure 
was also determined by the volumetric measurements and has a value of 
- 81 C. This value is lower than that obtained by the DSC measurements 
by about 22 C. The decrease in T^ is due to the different cooling rates 
(28 C/hr for the volumetric measurements and 1220 C/hr for the DSC 
measurements) and the different formation history used in each technique. 

The volumetric and light scattering transition data for 5P4E 
are shown in Figure IV- 22 and have behavior similar to that discussed 
above for Nl. 

From the above discussion, the glass transition determined from 
the light scattering technique is in general agreement with glass 
transition measurements obtained with other techniques. 

2 . Relatio nship of Glass Transition to EHD Conditions 

To obtain an impression of the relationship between glass transition 
for various lubricants and conditions in EHD contacts , Figure IV-30 has 
been prepared. Figure IV-30 contains the liquid-glass transitions for 
three lubricants as measured by the volumetric technique. The three 
lubricants selected represent a wide range of behavior with respect to 
glass transition. Above and to the left of each curve is the liquid 
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region for the material for low rate (long observation time) processes. 

To the right and below each curve is the glassy region. Based on 
limited measurements on three other hydrocarbon base oils the curve for 
N1 is an adequate first approximation for non-polymer containing 
hydrocarbon base oils. Also indicated in Figure IV-30 are the approximate 
ranges of pressure in the inlet and the Hertzian zones of an EHD contact. 
The approximate temperature ranges are also marked for an inlet temperature 
of about 20 C based on the IR temperature measurements reported elsewhere 
in this report. What this heuristic depiction shows is that glass 
transition most probably occurs in many practical applications of EHD 
contacts with many existing lubricants. The higher the contact pressure 
the more likely the transition to the glassy state, particularly with 
rolling and low slip ratio contacts. The glass transition and glassy 
state behavior is most likely to effect traction but may also effect 
film thickness if the transition occurs in the inlet zone, which can 

happen. The above clearly suggests the phenomena of glass transition 
occurs in EHD and its effects may be observable in EHD data already 
available. 

a. Traction Measurements 

Johnson and Cameron [27] measured the traction transmitted by 
an elastohydrodynamic film in a rolling contact with low sliding 
speeds and near isothermal conditions. From the traction measurements, 
the variation in apparent viscosity of their oil* with pressure and 
temperature is shown in Figure IV-31. The variation of the apparent 
viscosity with pressure shows a marked change above 10 Pas. The 

* 

Shell Vitra 76 similar to N1 but of higher viscosity. 
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authors [27] remarked that it is tempting to suppose that this change 
indicates a change in physical properties of the oil at high pressure. 
To test their remarks, the linear sections of Figure IV-31 were 
extrapolated to obtain an intersection. Pressures of about 0.57 and 
0.88 GPa were obtained at 30 and 70 C respectively. These intersection 
points fall near the glass transition lines of N1 fluid as shown in 
Figure IV 23 which is repeated with this data in Figure IV-32. The 
fact that the intersection points correspond to the measured glass 
transition points support the authors [27] remark, 
b. Side Thrust Measurements 

Johnson and Roberts [26] devised a rolling-contact experiment 
with a small amount of twist. Because of certain kinematic character- 
istics of their contact, no side thrust will be transmitted if the 
lubricant behavior is that of a liquid. If the material behaves like 
a solid, some side thrust will be observed. Figure IV-33 is a plot of 
the relative side thrust as a function of contact pressure. This 

transition in the behavior of their oil* is similar to that of specific 
heat or expansion coefficient of a material passing through a glass 
transition as discussed earlier. This observation raises the question 
of whether the behavior shown in Figure IV-33 is a glass transition. 

The onset of solid like behavior occurs at about 0.5 GPa and 23 C 
which is plotted on Figure IV-32. When the technique for determining 
the glass transition temperature for such a physical behavior was used, 
a midpoint transition pressure (see Figure IV-32) of about 0.7 GPa was 
obtained at 23 C. These values of pressure are seen to fall in the 
transition and glassy regions of Figure IV-32. This result agrees with 
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Johnson and Roberts' conclusion that a transition from viscous to 
solid behavior occurred above a contact pressure of about 0.5 GPa. 
c. Contact Temperatures and Pressures 

Fluid and ball surface temperature distributions in a sliding 
and rolling contact are shown in Section III of this report. Figure 
IV-34 is a plot of the ball surface temperature at the contact center 
versus maximum Hertz pressure at three sliding speeds. Figure IV-35 
shows the ball surface temperature versus sliding speed at different 
Hertz pressures. Because both figures have temperature and pressure 
plotted, we can include the glass transition lines obtained by the 
light-scattering and the volumetric measurement techniques. As shown 
in these figures, the lubricant near the ball surface is in the glassy 
state at a sliding speed of 0.2 m/s. At higher sliding speeds, the 
glassy state exists only at high pressures. This result was expected 
since a sliding contact generates frictional heat and thus the temperature 
levels in the contact are relatively high. 

The above disscussion suggests the existence of the glassy state 
in sliding EHD contacts at relatively low speeds or sufficiently high 
pressures. 


Ball surface temperature distributions for fluid N1 at different 

( u 2 ~ u x ) 

slide-to-roll ratio (Z = 2 ~z — rl are shown in Chapter III and are 

v^2 ' ^ 1 ^ ' 

repeated in Figure III-36. The speed ratio, Z, ranged from zero (pure 
rolling) to 2.0 (pure sliding) at a constant Hertzian pressure of 1 GPa 
and at constant rolling speed of 0.75 m/s. Based on the measured tempera- 
tures, assumed Hertzian pressure distributions, and the glass transition^ 
data, it can be estimated what portion of the fluid in contact with the 


70 



ball surface will be in the glassy or liquid states. The portions are' 
shown in the bar diagrams included in Figure IV-36. The lubricant near 
the ball surface is in the glassy state for all speed ratios except that 
of E = 2.0. At sufficiently low speed ratios, the glassy state occupies 
all the contact area. This result shows that the lubricant N1 is in the 
glassy state in a rolling/sliding contact and the glass can occupy all 
the contact area at low slide-to-roll speed ratios. 

G. Conclusions 

The glass transition temperatures of lubricants increase suffi- 
ciently with pressure to cause glass transition in EHD contacts with 
many typical lubricants. The glass transition behavior of lubricants 
of current interest varies considerably. Some, for example 5P4E, undergo 
glass transition at room temperature at about 0.14 GPa while bis-2- 
ethyl hexyl sebacate did not undergo a glass transition down to -3C 
even up to 1.1 GPa. 

It appears that EHD data already available in the literature 
are consistent with the glass transition observations reported here. 

Glass transition of lubricants can be expected to influence both the 
EHD traction and the film thickness depending on the location of the 
transition in the contact. 

The results of the preliminary study reported here clearly 
suggest that further research into the glass transition of lubricants, 
rheological behavior of lubricants in the glassy state and further 
verification of glass transition in EHD contacts are needed. Analytical 
modeling of the mechanics of EHD contacts with glassy lubricants 
should be pursued. Also the study of the relationship of chemical 
structure to glass transition may be useful in future lubricant 
formulation. 
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V. ON THE ANALYTICAL DETERMINATION OF PRESSURE 
DISTRIBUTION IN EHD LINE CONTACTS FROM FILM 
THICKNESS DATA 

A. Introduction 

The objective of the investigation reported in this section 
is the determination, by analytic methods, of the pressure distribution 
to which a film is subjected in an elastohydrodynamic process. The 
EHD-generated pressures in the film are the surface-tractions on the 
solid components of the system. In this study the solid components 
consist of either two elastic cylinders with parallel generators, or 
an elastic cylinder and an elastic plate. 

The interest inherent in the foregoing problem is many-fold. 
Firstly, the rheological properties of the film are functions of its 
local pressure and temperature. This dependence becomes more 
pronounced at the high pressures generated in the EHD-zone. Conse- 
quently the integrity of the film is tied to the understanding of this 
dependence. Secondly, if the intensity of the pressure field becomes 
large enough there is a distinct possibility that the lubricant’s physical 
state may change. Hence an accurate description of the thermodynamic 
variables is needed to assess the accuracy of the currently employed 
theories and hypotheses of the EHD process. Thirdly, the pressure 
distribution is necessary for the calculation of the stress fields to 
which the elastic components of the system are subjected. The deter- 
mination of the stress field could lead to a better understanding of the 
fatigue and fracture of bearing surfaces. Fourthly, since the deter- 
mination of the pressure field, as developed here, is independent of 
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surface forces 

elastohydrodynamic film thickness 

film thickness in undeformed state 

Fourier coefficient of pressure on half-plane 

half-plane solution corresponding to arbitrarily 
distributed surface load 

symmetric half-plane solution 

antisymmetric half-plane solution 

/-T 

set of real integers 
pressure 

Fourier coefficient of pressure 
reference pressure 
total applied load 
radius of cylinder 

radial location of arbitrary interior point z relative 
to points of application of surface forces 

residual 

antisymmetric cylinder solution 
symmetric cylinder solution 
sine integral 
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u,v 



U. 

1 


» 



V. 

l 


x,y 


X 

o 


z 

z k 

a 

a 

o 

3 

0 

Y 



A,y 

S 

a 

T ,T ,T 
xx xy yy 

<P,<P 

[ ] 


displacement in x- and y-directions respectively 
surface deformation of bodies 1 and 2 
coefficient of deformation 
cartesian coordinates 
reference length 

distance along x-axis measured from point of contact 
in the undeformed state 

complex variable x + iy 

points on surface of body at which forces act 
angular location of surface force 

angular half-width of pressure distribution on cylinder 
angular coordinate of arbitrary point on cylinder surface 
dimensionless angle, B/a Q 
normal approach of two solid bodies 

angular location of arbitrary interior point z relative to 
points of application of surface forces 

Lame elastic constants 

dimensionless coordinate x/a 

tangential surface traction 

stresses 

complex potentials 

stress tensor or displacement tensor fields associated 
with a particular solution 


Superscripts 


(a), (s) antisymmetric and symmetric solutions, respectively 

— — — complex conjugate function 

* differentiation with respect to the argument 
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the state of the lubricant, one may view the EHD-process as a reasonable 
method for examining the behavior of fluids under high pressures. 


B. Problem Statement 

ft 

As stated previously, the surface tractions , to which the 
elastic components of the system are subjected, are caused by 
the EHD generated pressure field in the thin film separating 
the solid components of the system. Since the region of the solid 
elements over which the EHD pressure is acting has dimensions which 
are very small by comparison to the dimensions of the solid, and 
furthermore the thickness of the fluid film is orders of magnitude 
less than its EHD-zone length, it is reasonable to assume that the 
variation of the EHD-pressure across the thickness of the film is 
negligible. 

Under these circumstances the surface tractions acting on 
either of the two solid bodies separated by the fluid film are 
identical. Figures V.l andV.2. Consequently, the problem under 
consideration, namely the determination of the pressure distribution 
in the EHD film, may be decomposed into three simpler problems. 

(i) The determination of the displacement field of the first 
solid body subjected to the EHD-generated pressure field over a region 
of its surface of the order of magnitude of the EHD-zone. 

(ii) The determination of the displacement field of the second 
solid body subjected to the same pressure field over the same region 
as the first solid body. 

*Surface tractions are the force distributions on the surface 
of the solid body. 
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(iii) The calculation of the film thickness as a function 
of position for a given pressure distribution. This task may be 
accomplished from the determination of the displacement fields of 
problems (i) and (ii) , the undeformed geometry of the elastic bodies, 
their normal approach and the total load transmitted across the film. 
The resulting relationships contain the film thickness and the EHD 
generated pressure intensity as dependent functions. The complete 
mathematical description of the EHD problem can be accomplished by 
deriving the appropriate field equations for the fluid film from the 
conservation equations of fluid mechanics. However, the coupling of 
the pressure field with the film thickness and consequently with the 
elastic displacements makes the resulting equation almost intractable. 
However, experimental techniques have been developed which make it 
possible to measure the fluid film thickness in an EHD process. Thus, 
assuming the availability of film thickness information, the complex 
complete mathematical problem may be bypassed in favor of the follow- 
ing simpler problem: the calculation of the pressure field in terms of 

the fluid film thickness. The resulting mathematical problem is much 
simpler than the complete one associated with the EHD process. 

In this investigation an EHD line contact is considered. Hence 
the elastic bodies bounding the film are either parallel cylinders or 
a cylinder and a plate. Section C contains the formulation of the 
boundary value problem of an elastic cylinder subjected to surface 
tractions over a small portion of its curved boundary, and the 
determination of the associated displacement field as a function of 
the EHD-generated pressure field. In Section D the deflection of a 
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thick plate subjected to normal loading is examined by approximating 
the plate by the elastic half-space and subjecting it to the same 
surface tractions employed in Section C. In Sections C and D it is 
assumed that the length of the cylinder is large enough so that the 
edge effects may be neglected and that the geometry of the bodies as 
well as the pressure distribution variations along the axial direction 
are sufficiently small to be ignored. With the adoption of these 
assumptions the elastic deformation problem reduces to one of plane 
elasticity and the complex function approach is employed. The coupling 
between the pressure distribution and the known film thickness is 
accomplished in Section E. 

The first phase of the problem is approached as a boundary value 
problem of the first kind (Sections C and D) . The result of the 
first phase is a relationship between the surface tractions and the 
film thickness. The measured film thickness is then employed in the 
second phase of the analysis to calculate the pressure distribution 
(Section E) . The numerical approach and its consequences are discussed 
in Section E. The influence of errors in the measurements of the film 
thickness on the pressure distribution and the method of minimizing their 
effect is also presented in Section F. Finally Section G deals with test 
cases for the evaluation of the whole analytic and numerical approach. 

C . The Displace m ent Solutions - The Cylinder 

In this section, the complex representation of the theory of plane 
elasticity is employed [64,65], for the purpose of formulating the solu- 
tions of an elastic cylinder under surface tractions. The surface 
tractions are assumed independent of the axial direction and the body 
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forces are ignored. Under these conditions, the stress tensor field 
and the associated displacement vector field may be represented in 
terms of two complex holonomic functions [66,67] by 


T xx + T yy = + ( c -l) 

T yy - T xx + 2iT xy = 2[z$"(z) + Y'(z)] (C.2) 

2(1(0 + iv) = $(z) “ zf'Cz) - Y (z ) (C.3) 

where ' implies differentiation with respect to the argument, denotes 

the complex conjugate function, and \,|j, are the Lam^ elastic constants. 

In the case of an elastic cylinder under surface forces F x + iFy, 
the functions $(z) and Y(z) are found to be 


§ ( 2 ) = — P — r 

) 2ni 3 5-z 2TTi J 


_ - a r z 

S ° 


*<*> = 2^i J ‘ - # ' (z > + - 


z z 


(0.4) 


(C.5) 


where 


a = a = $'(0) 



(C.6) 


If the surface forces consist of n concentrated forces applied at 
points of the circumference z^, the stress functions assume the form. 
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<P. 7) 


i(z) = 4^ L (F kx + if k y )[ 2 in(z k' z) ‘ z \ ] 

k=l 


n 


Y(z) = -t 1 [ (F kx - iF ky )to(z k' z > - < F kx + iF ky> ry < C ' 8 > 

k=l k 


Employing the above function for the particular case of a pair of equal 

and opposite forces, parallel to the real axis, acting at z l = r e ia , 

i (tt- of) , 

- re , one obtains 

z o 


$ (z) 




2 I 


(C.9) 


with 



the associated displacement field 

u = ■ nr ^ r- + c ° s 2e i - c ° s 2s 2 - (^) ~ cos 4 <c.n> 

^ O 

v = ‘ 4^ {(«{:) < e i + V - sln 29 i - sin 2e 2 ' (jSSr)f- cos “} ( c - 12 > 

o 


with p = F 2 = F ^ and (r^G^), signifying the location of an arbitrary 

interior point z of the cylinder relative to the points of application 
of the surface forces. The displacement field at any point on the surface 
may now be expressed in terms of the intensity of the load, p, the 
location of its application, (r Q ,o) and the arbitrary point on the surface 



u($;aO 


_ ?v+2 u J\ n cos %(o + ft ) 

2n(X-^)n, ^ sin%(a - f3) 


COSQ' cos 




(C.13) 


v(B! “ ) ” ' 2^0^' K“ + f) + X cos a sin g} (C.14) 

s 

This solution forms the basis for the subsequent development of the problem 
and is referred to as solution [B^]. The displacement field [ B ^ ] may be 
expanded by means of Green’s method, to yield the displacement field 
corresponding to an arbitrarily distributed surface loading of intens ity 
p(aO . 

In this work, however, the following procedure is adopted. Let [ S 2 2 J 
denote the solution corresponding to surface loading symmetric with respect 
to the x- and y-axes, and [S-^l denote the solution corresponding to 
loading antisymmetric with respect to the x-axis and symmetric with respect 
to the y-axis. Consequently the sought solution, [A], corresponding to 
the arbitrarily distributed surface loading is given by 

[A] = a^[S 12 ] + a 2 [S 22 ] (C.15) 


with a^ , a^ being constants and letters in brackets refer to the stress 
tensor or displacement vector fields associated with this particular 
solution. 

Solutions [S 12 ] and L s 22 J are subse< l uentl y generated by using the 
"basic" solution [B^] and employing Green's method with the appropriate 
choice of loading functions p (o) . 

1. The Symmetric Solution. 

Let two pairs of equal and opposite forces, parallel to the x-axis. 
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act at points z^, , z^ and z^ of the surface of the cylinder, with the 

points of application chosen in such a manner as to make the loading 
symmetric with respect to the x and y-axes. The displacement field 
associated with this loading is easily obtained from the basic solution 
[B 1 ], Figure (V-3a) 


u(3;a) = + 2 T t(\4u) pj^ntan 7r|a' + 3| + £ntan oi - 3 | + 20030^ os 3 j- (C.16) 


v(P;a) = ' (n+2o;) 


(C .17) 


The displacement field for distributed loading p(o0 over an arc (-or ,0^), 
which adheres to the symmetry imposed in this section, is 


A4.2U „ ^ r 1 1 

u(3) = tt ( )t-L J P( Q '){^ ntan + + Jfcntan * 3) 


O' 


+ 2 cosacosBtr do 1 


e}i 


O' 


v(s> = * ^0^0 I P<<*Kir + 2ar)rda 


(C. 18) 


(C . 19) 


Since the intensity p(ck') is symmetric with respect to the x-axis, which is 
the a = 0 ray, then p(o?) = pC-o). Hence it renders itself to the following 
represent at ion 


P(cO 


V /2n+l \ 

L P n COS l“ ™j 

~ O 


(G. 20) 
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which, in conjunction of Eqs . (C.18) and (C.19) and along with the 
assumption that the angles and 0 are small and the introduction [68,69] 
of 


q 


n 


2n+l 

2 


rr ; 



e 

S( 0 ,k 9 ) = J 

O 


s inkx , 
dx ; 


0 

c(o,ke> - J 

o 


1-coskx 

x 


dx 


(C. 21) 


yields the displacement field 


[S 


22 


]: 


■ I Pn U n (e; “o> 
n=o 

00 

v ») = £ p n v n (P!“ 0 ) 
n=o 


( C. 22) 


( C. 23) 


where 




a+2 “ )r o I q n“o . 

^<xV)n i ~T~2 Slnq n c ° SQ; o cosf3 
4 n o 


+ 


/O' v 

f 

n 


r _ a 1 

sinq n Xn(l+8) + j£n(l-3) + 2jfcn “M 


(C.24) 


+ sinqJ3[4n(l +3) - 4n(l - 3) + C(0,q n (l - 3)) 

- C(0,q n (l +B))] “ cosq n £[S(0,q n (l - 3)) +S(0,q n (l+ £))] !> 
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(C.25) 


V n (P ^o ) 


r a 
o o 

TT(\+|J.) 



+ n) 


sinq 3 
n 




2. The Antisymmetric Solution, [S.^] 

The loading in this case consists of two pairs of equal and opposite 

forces parallel to the x-axis, the a = 0 ray. The first pair acts at the 

points z, = r e iQ " and z„ = r e i ^ TT respectively, while the second pair 
r 1 o 2 o 

acts at = r e and z. = r and has its sense reversed with 

3 o 4o 

respect to the first pair. Thus the loading is antisymmetric with respect 
to the x-axis and symmetric with respect to the y-axis (Figure V-3b) . 

Taking advantage of the basic solution [B^, the displacement field 
associated with this type loading is found to be 


u(g;oO 


X+2u, 

2tT ( X4|Ji ) |JL 


p £n 


sin (crH3) 
sin(a-p) 


(C.26) 


v (3 ; cz) = - 


X 

tt (X+p)p 


P 


cos a sin(3 


(C . 2 7 ) 


Subsequently the displacement field for distributed loading over an arc 

(-a ,a ), with p^) = - p(-aO is simply 
o o 


u(3) 


X+2p 

2 TT (Xip, )|_L 



p(oO 


In 


sin(crfB) 

sin (cr*8) 


r^dcz 


(C .28) 


v(3) 


X 

TT (Xip)p 



p (o') cos CZ COS(3 


r dcz 
o 


(C.29) 


With the loading intensity p(cz) being antisymmetric with cz, its expansion 
is of the form 
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oo 


rmcy 


(or) = £ 


p sin 
n 


nrrcy 

a 


n=l 


(C. 30) 


Substitution of Eq. (C.30) into Eqs . (C.28) and (C.29), and adopting the 
notation introduced by Eqs. (C.21) along with assumption that and 3 are 
small, yields the results 


[S 


12 


]: 


u(0) 


l 


n=l 


p U (3; o' ) 
n n o 


(C.31) 


v(3) 


I 


p v (3; cy ) 

n n o 


where 


(C.32) 


(\+2jj,)r .a 

V 8; V = (^) f cos - -en(6+l)] + (cos q n 

+ cos q n B)[£n(3+l) - £n(3-l) - C(0,q n (j3-1)) - C (0 ,q^ (3+1 ) ) ] 


+ (sin q n £)[S(0,q n (l+S)) - S (O.q^l-B))]} 


(C .33) 


v (3 ;<* ) = 

n o 


Xr q 

o f n . \ 

rr(X+|jL)p, l 2 2 sin ^n/ 

q tck 


cosof o sin3 


(C.34) 


with q = nn. 
n 

3. The General Solution 

The displacement field corresponding to an arbitrary pressure 
distribution is obtained by employing the solutions [ ^22 anc * ^12^' 
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Its form is as follows 


U(P) 


= T {p< S) U< s >( B ;a o ) + P< 


a) 


n=0 



(C . 35 ) 


v«) = I K S)v n S)( S;« 0 ) +P^ a) V n a) (P^o>} < C - 36 > 

n=0 

where superscript (s) refers to the displacement field associated with 

the symmetric solution [S^], and superscript (a) refers to the displacement 

field associated with the antisymmetric solution [S^]. 

D. The Displacemen t Solutions - The Half-Space 

If one of the solid bodies that bounds the EHD zone from above or 
below is a thick elastic plate, whose dimensions are large by comparison 
to the EHD zone, then the elastic response of the plate can be represented 
reasonably by that of the elastic half space. If, furthermore, the other 
body is a sufficiently long elastic cylinder it can be assumed that the 
EHD-generated pressure is independent of the coordinate along the axis of 
the cylinder. Under these circumstances, the thick elastic plate can be 
substituted by the elastic half -plane. It is this case that is considered 
in this section and this is consistent with the formulation of the problem 
of section C. 

Again, as in the previous section, the displacement field may be 
expressed [66] in terms of two complex functions (j)(z) and i|^(z) by 

2p, (u + iv) = (~^) cp(z) - zcp T (z) - \|r (z) (D-l) 
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If the loading over the force boundary is 


T 

yy 

T 

xy 


= p(x) 
= ct(x) 


for x g [-a, a] 


it can be shown [66] that the stress functions <Kz) and iKz) are 
cp(z) = - 2 “ / (a + ip)j&n(t - z)dt 


H z ) = " ^ J (a + ip) + 2 tt J (cT “ ’ z > dt 


and the displacement field becomes 

4rrp,(u + iv) = - J* (a + ip) [(^^ L ) J^n(t-z) + Xn(t-z)J dt 
- (z-z) J* (a-ip) ~=r 

In this problem it is assumed that the tangential surface tractions 
are negligible (a = 0) , thus the displacement ea. (D.5) reduces to 

3 

4rrp,(u + iv) = - i J p(t) [(^^ a )^ n ( t "z) + j0n (t -z) Jdt 

-a ^ 

+ i(z-z) J p(t) 

-a 

Because any arbitrary loading function p(x) can be expressed as a 
sum of a symmetric and an antisymmetric component, the corresponding 
displacement field may be given as a combination of displacements 
corresponding to symmetric and antisymmetric solutions 


(D.2) 

(D-3) 

(D.4) 

(D.5) 

(D.6) 
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[H] = a[H 01 ] + b[H Q2 ] 


1. The Symmetric Solution, [H^] 

Assuming that p(x) = p(-x), p(x) may be expanded in the form. 
Figure V-4a 


P(t) = Y h cos 

Li m a 

n =1 


ra 


2n+l 
2 * 


n e J 


And the displacement becomes 


4 tt^ (v - iu ) 



^m «f ( cos £n(t-z)dt 

-a 

^cos S ^^^£n(t-z)dt 


+ 


2 i y 


oo 


Z h m I 3 

n=l 


^cos 


mrrt \ dt 
a / (t-z) 


Eq. (1).8), after some calculations and evaluating at y = 0, the free 
surface, yields 


4^ u = - Z ah m{< e i + ®2 )sin “m ' (*fr) t ( ®1 ' V sin q m 5 } 

. m 

n— 1 

00 

4 ^ v " - Z IT h m(x^ Ul ) { (sln q m )£n r l r 2 


(D.7) 


(D.8) 


CD. 9) 
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- (sin q m 5)[jtn ^ + C(0,q m (l+5)) - C(0,q m (l-?)j 


- (cos q m 5)[S(0,q m (l+5)) + S(0,q m (l-5))]} 


where 


9 1 +0 2 


0 for x > a 
TT for a > x > -a 
2tt for x < -a 


x 


§ = — , r, = a - x, r_ = a + x 


2n+l . ^ 

q = mrr = — x — tt as m section C. 


m 


2. The Antisymmetric Solution - [H Q1 ] 

The antisymmetric solution is characterized by Figure V-4b, 


p(x) = - p(-x) 


and consequently submits to the representation 


p(x) = 2, 


m s 


h sin 
m a 


m e J 


m=l 


Substitution of eq. (D.13) into eq. (D.6 ) yields 


(D . 10) 


(D . 1 1 ) 


(D.12) 


(D.13) 
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00 

4n^(v - iu) = - £ h m (sin M(t-z)dt 

1 3 

m=l 

00 

Cj i* (sin £n(t-z)dt 

m-1 * a V a ' 

+ 2iy I h m J_ ( Sln S! ?) (1% (D - l4 > 

1 ^ 
m=l 


Performing the indicated operations in eq. (D .14) and evaluating at y = 0, 
the displacement field at the free surface corresponding to the antisymmetric 
loading becomes 

00 

4TTM.U = - £ h m {^- (^)(cos <^5 - cos q^Ce^Bj)} < D ' 15) 

m 

m=l 


00 

4ttm.v = - £ h m ^ (^ Ul ) { (cos q m £ - cos qjtn ^ 
m=l 

+ cos q m § [C(0,q m (l+5)) - C(0 .c^Cl-g) ) ] 

- sin q m S [S(0,q m (l+§)) + S (0 .q^l-J) ) ]j- (D.16) 

3. The General Solution - [H] 

The displacement field corresponding to an arbitrary pressure 
distribution in x 6 [-a, a] is generated from the linear combinations of 
[H q2 ] and [H 01 ]. Thus 
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[H]: u(§) = £ 


(h (s) u (s) + h (a) u (a) ) 

~ m mm 


m m 


(D- 17) 


n=0 


v(s) = l 


( h (s) v ( s) + h (a V a) ) 

mm mm 


(D. 18) 


n=0 


( s ) 

where h is the pressure coefficient corresponding to the symmetric 
m 

loading, while h is the pressure coefficient associated with the anti- 
m 

symmetric loading, and 


U 


(s) a_ 

m 474J 


i (6 l + 9 2 )sln V - 


- 6 2 > sin V5j (D ' 19 > 


TT (s) a(A+2y) J, . 

V = - K- i r, . <(sinq JJlnr.r. 

m 2Trq y(A+y) | ms 1 2 


- (sin v 5) 


In + C(0 >qms (l + D) - C(0,q mg (l-O) 


J 


- (cosq ins OtS<0,q ins (l+C)) + S (0 .q^ <1-0 ) ] > 


(D. 20) 


U 


(a) _ 


m 


7-7-7 — r (cos q_ ^ - cos q ) (0 n - 0 O ) 

27rq ma (A+y) v \a M ma 7 1 2' 


„(a) a(X+ 2y) \ , _ r l 

V m = - 2rrq (X+y) \ (cos V 5 " cos V )iln 77 

ma ^ 2 


+ cosq C [C (0 , q ( 1+0 ) - C(0,q (1-0)] 

ma ma ma 


(D. 21) 


" sin ^ [S(0 ^a (1 ^ )) + S(0 ^a (1 - ? >^ 


(D. 22) 
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and 


q = mTT 
ms 


2n+l 


TT 


q = mTT - nTr 
ma 


n e J 


E. The Pressure Distributi on in the EHD Zone 

If the film thickness in the EHD contact is known from 
measurements, and the elastic properties and geometric characteristics 
of the solid bodies are known, it is possible to calculate the pressure 
distribution in the EHD zone. The film thickness may be expressed as 

h(x,p(x)) = h Q (x) - Y(P) + u 1 (x,p(x)) + u 2 (x,p(x)) (E.l) 

where h(x,p(x)) is the film thickness; h (x) is the film thickness in 

o 

the undeforraed state when the bodies are touching; Y (P) is the normal 
approach of the solid bodies; and u^ and are the displacements of the 
surface of the solid bodies at the EHD zone. The film thickness and 
the displacements are functions of position and of pressure distribution. 
The film thickness is known from measurements, and h Q (x) may be calculated 
from the geometry of the undeformed solid surfaces. In the previous 
two sections, expressions have been derived that give u^ and explicitly 
in terms of position and the Fourier coefficients of the pressure dis- 
tribution. It remains, then, to determine the normal approach Y(P) and 
the pressure coefficients. 

Specifically, since at any point A in the contact zone or its 
near neighborhood the film thickness is known, then 

h(x A ,p(x)) = h Q (x A ) - y(P) + U;L (x A ,p(x)) + u 2 (x A ,p(x)) (E.2) 
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where 


u l(x A , p (x)) 


N-l M f . f x 

J p (s) 0 (s) + I PmMm* 

L ‘ r n In L 


(E.3) 


n=0 


m=l 


u 2 <x ,p(x)) 


N-l M / \ / \ 

l + r c < 


(E.4) 


n=0 


m=l 


Csl) (s) 

IL , U, , are the symmetric and antisymmetric coefficients 

In ±m in ^m 

of the deformations of bodies 1 and 2 res p ectively , corresponding to a 
pressure distribution of the form 


N-l M 

V (s) 2n+l „ , v ( a ) j mTT /in 

p = ) p cos — = — ttx + ) p sin — x . (E.5) 

r u *n 2a L a 

n=0 m=l 

Further information is obtained by observing that the total pressure 
acting on either of the solid surfaces is equal to the applied load 
thus 


P 



p(x)dx 


(E. 6) 


If the film thickness is known at (N + M) points, there are (N + M + 1) 
equations available ((N + M) of the form of Equation (E.2), plus equation 

(’a) 

(E.6)) to determine the (N + M + 1) unknowns y(P) , p v , and p v . The 

n m 

pressure distribution in the EHD zone is then known from equation (E,5). 

It is apparent that if the EHD zone is bounded by two cylinders, 
both u^ and u 2 will be given by equation (C.35). If one of the boundary 
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surfaces is a cylinder and the other is a plate, is obtained from 
equation (C.35) and u^ from equation (D.18). The latter case is the 
one which will be dealt with in the following discussion; treatment of 
the case of two cylinders is analogous. 

Figure V-5 shows the conjunction of the cylinder and the plate. 

The pressure zone extends from x = -a to x = a, and 8 = -a to 8 = ot Q . 

Because a is very small, a = r a , and x = r 8. Since the location 
o o o o 

of the point x = 0 is the midpoint of the pressure distribion, 
it is necessary to specify the points at which the pressure goes to 
zero (i.e. the inlet and outlet points) in order to determined the 
location of x = 0 and the value of a. An ad hoc procedure for deter- 
mining these points is discussed in a later section. It should be noted 
that the origin of the x coordinate in this analysis in general is not 
the center of the dry Hertzian contact. 

F. Numerical Considerations 

As developed in the previous section, the problem of determining 
the pressure distribution in the EHD zone reduces to the problem of 
solving a set of algebraic equations for the Fourier coefficients of 
the pressure. Since equations (C.35) and (D.18) for the displacement 
of the solid surfaces are linear in these coefficients, the set of 
algebraic equations to be solved is linear. In order to avoid an 
ill-conditioned numerical system of equations, equation (E.2) is solved 
for the total displacement, 

u l(x A ,p) + u 2 (x A ,p) = h(x A ,p) - h Q (x A ) + y(P) A = 1,2, ... ,n (F.l) 
where the left-hand side is expressed in terms of the Fourier coefficients 


93 



P iS i = 1,2,... t m by equations (E.3) and (E.4). A trial value of y(P) 
is assumed, and when m = n the system (F. 1) may be solved for p^. 

Equation (E. 6) is then evaluated to see if the total load condition 
is satisfied. If not, a different value for y is selected and the process 
repeated. Since the system is linear, two iterations are required to 

find the true value of y. Because only the right-hand side of equation 
(F. 1) is affected by changing y, very little extra work is incurred by 
this iteration scheme. If the film thickness measurements were known 
to be exact, the film thickness at m points would be sufficient to 
determine the Fourier coefficients. However, due to unavoidable 
inaccuracies in the measurements, certain modifications to the solution 
of the equations must be made. 

Since the measurements at each point in the EHD zone are in error, 
the total error increases as the number of points used, n, increases. 

At the same time, a sufficient number of points must be used to 
adequately characterize the shape of the film thickness profile. 
Furthermore, the number m of Fourier coefficients to be found must be 
large enough that the truncated Fourier series, equation (E.5) (with 
m = N + M) may accurately represent the true pressure. To alleviate 
this problem, n and m are chosen, with n > m, and the following process 
used to "smooth out" the errors in the film thickness. The algebraic 
equations take the form 

a^. p. = b_£ i = 1,2 ,. . .n, j— l,2,...,m (F.2) 

where p^ are the unknown coefficients. Equation (F.2) represents an 
overdetermined system (that is, there are more equations than unknowns) 
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for which, in general, there is no exact solution. Accordingly, the 
residual, r^, is defined as 


a ij p i - b ± = r ± 


(F.3) 


The unknowns are then determined so as to minimize the sum of the squares 
of the residuals, r^r^. This results in a system of m linear equations [70], 
The minimization of the residuals has the effect of smoothing the data 
while coming as close as possible to solving equations (F.2), 
recognizing the fact that the data is a best— estimate to the 
true film thickness. The amount of smoothing depends on the degree to 
which n is greater than m. If the Fourier series for the pressure is to 
converge to the true pressure, the coefficients must eventually go to 
zero as m is increased. Once n is selected, the amount of smoothing 


required is determined by increasing n until the Fourier coefficients 
begin to increase, causing the series to diverge. As a check on the 
acceptability of the result the magnitudes of the residuals may be 
compared to the expected error in film thickness measurements. 

It should be noted that the calculations for the pressure 
are based on the displacement of the surfaces from their undeformed 
states, rather than on the film thickness Itself. From equation (E.l), 
assuming h^Cx) and Y(P) may be found precisely, the magnitude of the 
error in total displacement, u^ + is the same as that for the film 
thickness. Consequently, near the center of the contact where the film 
thickness is small and the surface deformation is large, a small relative 
error in film thickness results in a smaller relative error in the 
displacement. Therefore, the pressure calculated in the part of the 
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contact where deformations are large is expected to be quite accurate. 

On the other hand, in the inlet region where the pressure is building 
up gradually from zero, the film thickness is large due to the initial 
surface geometry while the displacements are small because of the lower 
pressure. Hence, a small relative error in film thickness in the 

inlet region may result in a large relative error in the displacement, 
so that the pressure calculated in this region is somewhat less reliable. 
Thus it is often helpful to introduce additional information in this area 
by requiring that the pressure gradient at x = -a vanish. This is a 
consequence of the fact that the pressure is generated by the hydro- 
dynamic action of the fluid, and has the effect of introducing an 
additional equation into the set to be solved. This procedure tends to 
offset somewhat the increased effect of errors in the inlet region. 

Because of the method by which pressures are calculated, an 
error in film thickness which is constant in magnitude throughout 
the entire EHD zone has no effect on the computations. Such an error 
has the character of a rigid-body translation of one of the solid 
surfaces , and would therefore be absorbed in the normal approach 
term (y) of equation (E.l), leaving the displacement and hence the 
pressure unaffected. Consequently the range of errors encountered, 
rather than their magnitude, determines the accuracy of the calculated 
pressure distribution. 

G. Test Cases and Evaluation 

In order to test the procedures described in the previous 
section, and to evaluate the effects of certain types of errors in 
the input data, a number of test cases were run. In the first case, 
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the following pressure distribution was employed 


p(x) 



(G.l) 


where x = 0 is the contact point of the undeformed surfaces. Closed-form 
solutions for the displacements of the cylinder and the plate for the 
pressure distribution were found from equations (C.13) and (D. 6). 

These displacements were then used as film thickness data to obtain a 
series expression for the pressure of the form of equation (E.5), which 
could be evaluated by comparison with equation (G.l). The calculated 
pressure curve which resulted is shown in Figure V-6, which represents 
equation ( E .5) with 10 terms, and a = 1.3 x^. The major disagreement 
with equation (G.l) occurs near £_= + x q where the pressure is discon- 
tinuous; this is to be expected due to the nature of the series, and is 
not a serious drawback, since pressure discontinuities are not expected 
in a film. The selection of a = 1.3 x q results in small negative 
pressures near the inlet and outlet. This will be discussed later in 
this section in more detail. 

As a second case, the film thickness was taken as zero over the 
region - x q < x < + x q ; this is the case of two dry bodies in contact 
over the region, for which the pressure distribution should approximate 
the Hertzian semi-elliptical distribution. The results are shown in 
Figure V-7 for nine terms in the series of equation (E.5), with a and 
x = o chosen so that the pressure is zero at x = - 1.33 x and x = 1.1 x . 
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The agreement is good except near x = + x q , where there is a discontinuity 
in the pressure gradient. In the presence of such a discontinuity, 
inaccuracy in the series expansion is to be expected, since equation 
(E.5) implies continuous derivatives. 

A third test case was examined using a pressure distribution 
expressed by a third-degree polynomial, shown in Figure V-8. As in the 
first case, a closed-form expression for the displacement of the surfaces 
was obtainable, and the resulting film thickness profile is also shown 
in Figure V-8. The values of the total load, elastic constants of 
the surfaces, cylinder dimensions, and minimum film thickness were 
selected so as to approximate those in an actual EHD contact [71]. 

The calculated film thickness was used to obtain a series 
solution for the pressure as described in the preceding section. When 
the inlet and outlet points (x = + a) were specified correctly, the 
series expression for the pressure showed negligible deviation from the 
polynomial expression throughout the contact. Because the pressure 
distribution of Figure V-8 in some ways resembles that expected in an 
EHD contact, this test case was examined further in order to study the 
effects of certain types of errors. It was estimated that, using the 
optical interference method of measuring film thickness, the maximum 
error in film thickness at a given location is .025 ym (1 yin.). 
Accordingly, three sets of random errors in the range (-.025, .025) ym 
were introduced into the film thickness data. This was a fairly severe 
test, because in actuality the errors would most likely be normally 
distributed, over the range (-0.025, 0.025 ym) . The pressure distribu- 
tions calculated from these three error sets are shown in Figure V-9. 
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Sets 1 and 3 show the largest errors on the side corresponding to the 
inlet of an EHD contact, as would be expected from the discussion of 
the previous section. The general shape of the pressure distribution is 
preserved in all three cases shown in the figure. 

The effects of errors in location of the beginning and end of 
the pressure were also studied. The results of varying the point of 
zero pressure at the inlet of the EHD zone are shown in Figure V-10. 

When the inlet point is located too far in, a hump in the calculated 
pressure appears, and convergence of the series becomes somewhat slower. 
The presence of these conditions indicates that the inlet point should 
be moved farther out. When the inlet point was located too far out, 
the calculated pressure came in at zero up to the true inlet, and 
blended into the exact distribution almost immediately, resulting in 
negligible error. When the change in pressure near the inlet is more 
abrupt, as in Figure V-6, negative pressures may result if the inlet 
point is set too far out. Since the location of the point of zero 
pressure is not known ahead of time, these conclusions are helpful in 
locating the inlet point when calculating the pressure distribution. 
Similar conclusions may be drawn concerning the location of the outlet 
point; however, the determination of the location of this point in an 
actual EHD contact is generally less difficult, since it is known a 
priori that the pressure goes to zero in the vicinity of the exit 
constriction of the film. 

In addition, small (3%) errors were introduced in the elastic 
constants, the total load, and the radius of the cylinder. The cal- 
culated pressure distribution was found to be relatively insensitive 
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(- 1%) to these errors in the elastic constants and the load, and 
somewhat more sensitive to errors in the radius. This is not a serious 
drawback, since the radius may be measured to a high degree of precision. 

At present, accurate data for the film thickness in an actual 
EHD line contact are not available to us. Wymer and Cameron [71] show 
measured line contact film thickness profiles, but from their published 
figures it is possible to determine the film thickness only to within 
2-3 |i in. at a given location in the EHD zone. Nevertheless, a pressure 
distribution was calculated based on their measurements, and the results 
shown in Figure V-ll where it is compared with the Hertzian distribution. 
The figure also shows a comparison between the measured film thickness 
used as input for the calculations, and the film thickness which would 
result from the calculated pressure distribution. The difference between 
the two is the residual, which is quite large near the exit region, as 
shown in the figure. In this region, the magnitude of the slope of the 
film thickness profile is quite large, which increases the difficulty 
in reading the points from the figures. Because of the large residuals 
near the exit and the inaccuracy of the input data used, the pressure 
distribution of Figure V-ll should be viewed as a rather rough approxi- 
mation. In spite of the lack of precise film thickness data for a line 
contact, accurate data is available for the EHD point contact. Most of 
the foregoing numerical methods and conclusions concerning probable 
errors are expected to carry over into this case. 

H. Summary of Computational Method 

If the film thickness is known from measurements at n points 
A = l,2,...,n, in the EHD zone, the set of n equations (F. 1) is used 
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to find the Fourier coefficients of pressure as follows. The film 
thickness in the undeformed state is known at each point x^, and an 
initial estimate for the normal approach y(P) is made. The right- 
hand side of equations (F.l) is then known. The displacements u^(x^,p) 
and u„(x & ,p) are expressed in terms of the m unknown Fourier coeffici- 
ents. The number of Fourier terms used (m) is arbitrary as long m < n. 

For the case of two cylinders, both u^ and u^ are given by equation (C.35); 
for the case of a cylinder and a plate, u^ is given by equation (C.35) 
and u 2 by equation (D.18) . Equation (F.l) is then a set of n equations 
in m unknowns. An additional equation is obtained from the requirement 
that the pressure gradient be zero at the inlet (4^| _ = 0 from equation 

(E.5)). Because the system is overdetermined (m < n + 1) residuals are 
defined as in Equation F.3. The condition that the sum of the squares 
of the residuals be a minimum is applied, which results in a system of m 
equations in m unknowns. This system is solved for the m Fourier coeffi- 
cients, for use in equation (E.5) for the pressure distribution. The 
total load P is calculated from equation (E.6), and the result compared 
with the known value of the load. If they are not the same, a new estimate 
is made for y(P) , and the solution process is repeated. Once equation 
(E.6) is satisfied, the Fourier coefficients of pressure are known, and 
the pressure at any point in the EHD contact is determined from equation 
(E.5). 

The calculations require knowledge of the film thickness at n 
locations in the EHD zone, the total load on the system, the elastic 
constants of the solid bodies, and the radii of the cylinders. In 
addition, the number of sine and cosine terms to be used in equation 
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(E.5) and the locations of the inlet and outlet points must be 
specified. On the following pages is a listing of the Fortran 
computer program for calculating the pressure distribution by the 
above method. 
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77/09/11. 10.13.50. 

GEORGIA '''ECH CYQE^ 1 74. NO? l.t-430 

TTSRn MTTMBEPs ME650AC 
PASSWORD 

mum 

'TERMINAL: 77# ' r * r Y 

RECOUP /^vc'i'pi. BATCH# 50000 
$^FL# 50000. 

/GE^MLSDP 
/LNH# F=MLSDR 

PROGRAM EHD n ( DA' 7 '# I N n UT # OUT 13 TIT# TAPE 5- DAT# TAP E4= INPUT# TAPE 6- OUTPUT) 

C nn 0 C n AM TO COMPUTE n RES SURE DISTRIBUTION FOR EHD LINE CONTACT 

C FROM FILM THICKNESS DATA. 

DIMENSION H ( 30 ) # BETA 1 < 30 > # BETA2 ( 30 ) # 0 ( 30 > # Al ( 30# 30 > #XT( 30 > 
DIMENSION A2(30#30)#AS(30#30)#Y(30)#UC30)#PC51 ) # YIN < 5 !) # XTI N < 5 1 ) 
DIMENS I ON I "I UO' 1 * C 30 ) # D ( 30 ) # SL C 30# 30 ) # UI N < 30 ) 

DIMENSION UC ( 30 > # HC ( 30 ) # RES C 30 ) 

REAL K1#K2 

C PRELIMINARY DATA: 

C SET I FLAG=0 FOR T«0 CYLINDERS# =1 FOR CYLINDER + PLATE 

C GAM: INITIAL CHESS FOR NORMAL APPROACH OF BODIES CM) 

C F: TOTAL LOAD CN) 

C M: NO. OF COS TERMS IN SERIES FOR PRESSURE 

C N: NO. OF SIN TERMS IN SERIES FOR PRESSURE 

C LI: NO. OF FILM THICKNESS DATA POINTS TO BE INPUT 

C sp; DESIRED SPACING BETWEEN hqintS AT VHICH PRESSURE 

C IS TO BE PRINTED OUT CM) 

C Y1#Y2: INLET# 0TTTLE' r ^OINTS# MEASURED FROM POINT OF CONTACT 

C IN TRIDEFORMED «TATE CM) 
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FEADC5**) I FLAG* GAM * F 
READC5* *>M*M*L1 * sp 
P EAD( 5 ,*)Yl,Y2 
p I “4 • *ATAN ( 1 . ) 

A=ABS( (XI -V2 ) /2 . ) 

0FST=-(YI +X2>/2. 

NN= 0 
MM=M- 1 
L 1 l -L 1 + 1 
LL=M+N 

FILM T HICKNESS DATA: 

X( I ) : DISTANCE OF I -TH DATA POINT FPOM POINT OF CONTACT IN 
TTNDEFOPMED STATE CM) 

H ( I ) ! FILM THICKNESS AT XCI) CM) 

P EAD<5**> (XCI )jH(I )* 1 = 1*L1 ) 

DO 13 I-1*L1 
YT< I )=:'( I ) + OFST 
13 CONTINUE 

DATA ON BODY 1 - CYLINDER: 

YM1: MODULUS OF ELASTICITY (N/M**2) 

PP 1 : °0I S.S0N ' S P ATI 0 
nis RADIUS CM) 

p EADC5**>mi*P p l*Rl 

K 1 3 2 . * ( ! . - pp 1 * pp 1 )/ ( D I *YM 1 ) 

ALP 1=A/P I 
DO 5 I - 1 * L 1 
BE^Al ( I )=XT( I )/ p l 
5 CONTINUE 

CALL CYL ( ALP 1 j 3ETA1 *N*M*L1*A1) 

I F ( I FLA G • EQ • 0 ) GO TO 10 
C DATA ON B0D v 2: USE IF I FLAG= 1 

PEADC 5j * ) v M2* pp 2 
K2= 2 • * ( 1 .- p R2* pp 2)/C p I*YM2> 

CALL PLATE<A*XT.rN*M*LI * A2) 

”1^0 .0 
GO ^O 17 
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C DATA ON BODY 2s USE IF IFLAG=0 

10 READC5, *>YM2,PR2,R2 

K2=2 • * < I • -PR2*PR2 ) /(PI *YM2 ) 

ALP2=A/R2 
DO 20 1=1 /LI 
BETA2CI )=X?CI )/R2 
20 CONTINUE 

CALL CYL C ALP2, BETA2, N, M, L l , A2 ) 

RIN=1 ./P2 

COMBINE MATRICES FOP BODIES 1 AND 2. SLCI,J> IS THE MATRIX OF 
COEFFICIENTS OF THE UNKNOWN FOURIER COEFF. OF PRESSURE. 

17 DO 25 I = 1 , L l 
DO 25 J=1,LL 

SLC I , J ) = C A1 (I> J)*K1+A2CI, J)*X2)*A/pI 
25 CONTINUE 

INCLUDE CONDITION THAT PRESSURE GRADIENT VANISH AT INLET 
DO 70 I = LN 

SLCLl 1, I ) = ( <- 1 • )**I*I >*( 1 .E-10) 

70 CONTINUE 
DO 71 1=1, M 

SLCLl l,N + I >=C C- 1 . )**CI-1 >*C2.*I-1 . )/2. )*1 .E-IO 

71 CONTINUE 

FROM SL, FORM MATRIX CORRESPONDING TO MINIMIZATION OF THE NORM 
OF THE RESIDUAL. 

DO 55 1 = 1, LI- 
DO 55 J= 1 , LL 
ASCI, J) =0.0 
DO 55 K= 1 , L 1 1 

ASCI, J)=ASCI, J)+SLCK, I )*SLCK, J) 

55 CONTINUE 

3EGIN PROCEDURE FOR SOLUING THE SET OF LINEAR EQUATIONS FOR THE 
UNKNOWN FOURIER COEFF. 

DO 1 1=1, LL 
I PIVOT C I )=I 
R0 rT MAX=0 • 0 
DO 2 J= 1 , LL 
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2 R 0 ^TM AX- AM AX 1 ( RQ T 7MAX.* ABS CASCI/J) ) ) 

1 D(I)=sRCWMAX 

NM 1 -LL- 1 
DO 3 K=»1#NM1 
J=K 

KP 1 =K+ 1 
IpalPI W0T<K) 

COLMAX- ABS (AS(IPjK) ) /D ( I P > 

DO A I = K n 1 » LL 
IP=IPIU0T< I ) 

At : tIKO^=ABSCAS( IPj K) 5/DCIP) 

IFCAUIKOU.LE.COLMAX) GO TO 4 

COLMAX=AUlKOU 

«J= I 

4 CONTINUE 

I °K= IPI POT ( J ) 

I P I V OT ( J ) = I P I T , T 0 ' T> C K ) 

I PI^OT (K)“ I PK 
DO 3 I»Kt>I,LL 
IP=IPIVOT< I ) 

AS ( I Pj K ) = AS ( I P* X) /AS < IPKj K) 

RATI 0=-*ASCIPj>K) 

DO 3 J=K©l^LL 

3 ASCI^jiJJ^RA^I 0*A$ ( I PTO J)+AS(IP.»J) 

C FORM RIGHT-HAND SIDE OF THE EONS* 

90 DO 45 I = 1 .» H 

UINCI )=.5*X(I >*XCI )*C 1 ./R1+RIN)-H< I ) - GAM 
45 CONTINUE 

FROM UIN AND SL> FORM RHS OF EONS* CORRESPONDING TO MINIMI 
ZATION OF THE NORM OF THE RESIDUAL. CONTINUE SOLN. OF 
LINEAR EONS. WITH THIS AS RHS. 

T.TINCL1 1 >=0.0 
DO 36 I - 1 !/ LL 
TT(I 5=0.0 
DO 36 . J- 1 *L 1 1 
U(I ) = tt<I >+UIN(J)*SL< I ) 
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36 CONTINUE 

I T3= IP J VOT ( i ) 
qci )=ucip) 

DO 9 K=2j LL 

IP=IPIttot(K) 

KMI-K-I 
S TI M=0 *0 
DO 6 JaUKMl 

6 STJM-AS(IP* J)*Q<J)+SVM 

9 0 ( K) =V ( I P ) - 5 TIM 

0 (LL )=0 (LL ) /AS < I°j LL) 

K=LL 

DO 7 NP I MK-2 > LL 
KP1 =K 
K=K- ! 

I P* IPIVOTCK) 

?TTM = 0.0 
DO 8 J=KPWLL 
8 ST1M=AS< IP* J)*Q( J) + STTM 
C OCK)! DESIRED FOURIER COEFFICIENTS 

7 Q(K5*(QCK)-STJM>/AS( IP,K) 

C IF NECESSARY* MODIFY VALUE OF GAM TO SATISFY TOTAL LOAD CON- 

C DITI ON* AND RE-SOLVE THE EQUATIONS VI TH HEW RHS* 

FT°I AL=0 • 0 
DO 34 1 3 1 * M 

T=4.*A*(-1 .)**(I-l)/((2.*I-l . >*PI ) 

34 FTRIAL=FTRI AL+T*Q<N+I ) 

FO=FTPI AL-F 
NN=MN+ 1 

I FCNN • EO • 3 ) GO TO 100 
I F ( NN • EQ • 2 ) GO TO 110 
F I = F0 
GAM 1 = GAM 
GAM= GAM/2. 

GO TO 90 



no GAM 2= GAM 1 
GAM 1 a GAM 
F2= Fl 
F 1 = FO 

GAM- GAM 1 - ( (GAMI -GAM2)*F1 )/<Fl-F2> 

GO TO 90 

COMPTTTE pRESGTtre AT LOCATIONS DESIRED FOR OUTPUT. 
100 NP= I NT CABS CXI /SP ) ) +INTC A3SCX2/SP) )+3 
VTINC 1 ) = AM INI (XI ,X2) 

XT I N ( 2 ) = I NT ( XT IN(1)/SP)*SP 

N Pl- N p-l 

DO AO I = 3.» NP 1 
XTINC I )*XTINC I- 1 ) +SP 
AO CONTINUE 

XTI N (NP ) =AMAX1 <X1,X2) 

DO A2 L S LN? 

XIN(L)=XTIN(L) +OFST 

PCL)=0.0 

DO 32 I = LN 

P (L )=P (L) +0. ( I >*SIN< I*PI*XIN(L)/A) 

32 CONTINUE 

DO 3 5 JJ=nM 
I=JJ- 1 

P(L)=p(L)+QCN+JJ)*C0S( C2.+I+1 . ) *PI *XIN (L)/ ( 2. *A) ) 
35 CONTINTTE 
A2 CONTINTTE 

COMPUTE RESIDUALS AT EACH DATA POINT LOCATION. 

TP= 0 . 0 

do 9i i = mi 

TJCCI > = 0.0 

DO 92 J=1,LL 

UCCI )=UC(I )+SL< Ij J)*Q( J) 

92 CONTINUE 

HC < I >= • 5*X ( I )*X< I )*< 1 ./P.1 +RIN)-UC ( I ) - GAM 
RES (I)=H(I)-HC(I) 

TR=TR+RE5 < I )*RES ( I > 
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91 CONTINUE 
C OUTPUT: 

C FTRIAL: CALCULATED TOTAL LOAD CN) - USE FOR COMPARISON AS A TEST 

OF ACCURACY OF THE CALCULATIONS. 

GAM: NORMAL APPROACH VALUE ARRIVED AT CM) 

FOURIER COEFF.: COEFFICIENTS OF SIN TERMS LISTED FIRST 
RESIDUALS: LISTED IN SAME ORDER AS INPUT POINTS CM). 

LAST VALUE PRINTED IS SUM OF THEIR SQUARES- 
X: DISTANCE FROM CONTACT POINT IN UNDEFORMED STATE 
VP I TEC 6, * ) FTR I AL, GAM 
VRI TEC 6, 1 50) 

150 FORMAT C 5X, "FOUR I ER COEFFICIENTS:”) 

VRITEC6, *) CQCI ), 1=1 AL) 

VRI TEC 6, 250) 

2 50 FORMAT C5X, "RES I DUALS:") 

VRI TEC 6**) CRESC I ), 1=1, LI ) 

VRITEC6,*)TR 
VRI TEC 6, 550) 

VRITEC6, 450)CXTINCI ),PCI),I=1,NP) 

550 FORMATC 10X,"X CM ) ", l OX, "PRESSURE CN/M**2)"/) 

450 FORMAT C 9X, F 1 0 • 5, 1 OX, El 1 *4) 

VRI TEC 6, 650) 

IF USER WANTS CALCULATED FILM THICKNESS AT EACH INPUT POINT, 

ENTER l. OTHERWISE ENTER 0. 

650 F0RMATC5X, "ENTER 1 FOR FILM THICKNESS") 

READC4, * )NP 
IFCNP.NE. 1 )G0 TO 93 
VRI TEC 6, 750) 

750 FORMAT C l 2X, "X CM) ”,9X,"FILM THICKNESS CM)"/) 

VRI TEC 6, 450) CXC I ) ,HC C I ) , I = 1 , L 1 ) 

93 VRI TEC 6, 8 50) 

C IF USER WANTS CALCULATED FILM THICKNESS AT OTHER POINTS IN OR 

C NEAR CONTACT, ENTEP 1. OTHERWISE, ENTER 0. 

850 FORMATC 5X, "ENTER 1 FOR ADD *L FILM THICKNESS") 

READ C 4, * )NR 
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IF<NP.NE. 1 ) STOP 

C ENTER LOCATION OF DESIPED FILM THICKNESS CALCULATION. TO STOP PROGR- 

AM.* 

C ENTEP VALUE > 20. 

94 READ ( 4* * )XX 
IFCXX.GT.20. ) STOP 
XT ( 1 )=XX+OFST 
BETA1 (1 )=X?< 1 >/*M 

CALL CYL< ALP 1, BETA 1,N,M, 1*A1) 

IFC I FLAG. EQ .0 ) GO TO 95 
CALL PLATE(A,XT,N,M, 1*A2) 

GO TO 96 

95 BETA2 ( 1 ) =XT ( 1 )/R2 

CALL CYL<ALP2,3ETA2,N,M, bA2) 

96 DO 97 1*1, LL 

SL< 1 , I ) = <A1 (1,1 ) *K1 +A2 ( 1, I )*K2)*A/PI 

97 CONTINUE 
UA*0 *0 

DO 98 1*1, LL 
UA-TJA+SLC 1, I >*Q(I ) 

98 CONTINUE 

HH* . 5*XX*XX* ( 1 • /Ri +RIN) -UA-GAM 

VR.ITEC6, 450 )XX, HH 

GO TO 94 

STOP 

END 

SUBROUT I NE CYL < ALP, BETA, N, MT, L 1 , A ) 

C SUBROUTINE TO FILL MATRIX CORRESPONDING TO CYLINDER DISPLACEMENTS 

DIMENSION A(30,30), BETA ( 30 ) 

REAL M 

PI = 4.*ATANU • ) 

MM*MT-1 
LL=MT+N 
DO 10 L= 1 ,Ll 
NN=0 
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I FCABSC 1 .-ABS (BETA CL) /ALP) > #LT. I *E-10)NN=1 
C CALCULATE ELEMENTS CORRESPONDING TO ANTISYMMETRIC COMPONENTS 
DO 12 1=1, N 
COEFl= I *PI /ALP 
XL=COEF 1 * (BETA ( L ) +ALP ) 

CALL SICOCXL, SI ,C01 ) 

XL=COEFl * (BETA (L) -ALP) 

CALL SI CO (XL, S2, C02 ) 

IFCNN.NE* 1 )GO TO 16 

A1=0.0 

GO TO 15 

16 Al=( (-1 . )**I-C0SCC0EF1*BETACL) ) ) *ALQG(ABS( (BETACL) 

1 +ALP)/ (BETA(L) -ALP) ) ) 

15 Bl=SIN(COEFl*BETA(L) )*( S2-S1 ) 

Cl=COS(COEFl *3ETA(L ) )*(C01-C02) 

ACL, I )=( 1 ./I )*(A1+B1 +C1 ) 

12 CONTINUE 

C CALCULATE ELEMENTS CORRESPONDING TO SYMMETRIC COMPONENTS 
DO 14 J J= 1 ,MT 
I=JJ-1 

M=(2.*I + 1 0/2. 

C0EF2=M*°I /ALP 

XL - C 0 EF2 * ( B ETA ( L ) + ALP ) 

CALL SICOCXL, SI ,C01 ) 

XL=C0EF2* (ALP -BETA (L> ) 

CALL SICOCXL, S2,C02) 

A2= ( 2 • *M* (PI ) **2/ CM**2*PI **2+ALP**2 ) ) * 

1 SIN(M*PI ) *C0 S (ALP )*COS( BETACL) ) 

IFCNN.NE. 1 )GO TO 20 

B2=SIN(M*PI )*(2.*AL0GC2.*ALP)-AL0G(4. ) ) 

C2=0.0 
GO TO 25 

20 B2=SIN(M*PI )*(ALOGCABS( (ALP+BETACL) )*(ALP- 
1BETACL) ) ) )-AL0GC4. ) ) 

C2=SIN(C0EF2*BETA(L) )*(ALOG(ABS( (ALP+BETACL))/ 

1 (AL°-BETA(L) ) ) ) > 



25 D2=SIN(C0EF2*BETACL) )*<C02-C01 ) 

E2= ( - 1 • )*C0S(C0EF2*BETA(L) >*< S2 + SI > 

A(L*N+JJ) =A2+ 1 • /M* (B2+C2 + D2 + E2 ) 

14 CONTINUE 
10 CONTINUE 

RETURN 

END 

SUBRO UT I NE PLATE ( A* X, N > MT* L \> AN ) 

C SUBROUTINE TO FILL MATRIX CORRESPONDING TO PLATE DI SPLACEMENTS 

DIMENSION AN<30,30),X(30> 

REAL M 

P 1=4 • * AT AN < 1 . > 

MM=MT-i 
LL=MT+N 
DO 10 L=1,L1 
NN= 0 

I FCABS C 1 • -ABS (XCL ) /A) ) *LT • 1 • E- 1 0 )NN= 1 
C CALCULATE ELEMENTS CORRESPONDING TO ANTISYMMETRIC COMPONENTS. 

DO 12 I = 1 > N 
C0EF1=I*PI/A 
XL= C 0 EF 1 * C A+X ( L ) ) 

CALL SI CO (XL* S 1 * CO 1 ) 

XL=C0EF1*(A-X<L) ) 

CALL SIC0CXL,S2.,C02) 

I FCNN.NE* 1 ) GO TO 16 

A1=0.0 

GO TO 15 

16 Al=( Crl • >**I-C0SCC0EF1*X(L> ) ) *ALOG( ABS ( ( A+X CL ) ) / 

1 (A-XCL) >) ) 

15 B1 = ( - 1 . )*SINCC0EF1*XCL> )*<S1+S2> 

C 1 = C0S (C0EF1 *X(L) )*(C01-C02) 

ANCL^ I)=(1./I)*(A1+B1+C1) 

12 CONTINUE 

C CALCULATE ELEMENTS CORRESPONDING TO SYMMETRIC COMPONENTS- 
DO 14 JJ=1,MT 
I=JJ~1 
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M»<2.*I+1 . )/2. 

C0EF2=M*PI/A 

S=SINCC0EF2*XCL>> 

XL-C0EF2* (A+XCL ) ) 

CALL SIC0CXL,S1,C01 ) 

XL=C0EF2*<A-X<L) ) 

CALL SI CO (XL, S2, C02 ) 

IFCNN.NE. 1 )GO TO 20 
A2-AL0GC2. )*C2.*SIN(M*PI ) ) 

B2=0.0 

GO TO 25 

20 A2-AL0GCABS ( < A-X(L) )/A) ) * < SIN <M*PI ) -S > 
B2“AL0G(ABS ( CA+XCL) )/A) >*< SIN(M*PI ) + S) 

25 C2= (- 1 .)*S*(C01-C02) 

D2= C - 1 . >*C0S(C0EF2*X(L) >*(S1 +S2) 

AN(L, N+ JJ ) = ( 1 • /M ) * ( A2+B2+C2+D2 ) 

14 CONTINUE 
10 CONTINUE 
RETURN 
END 

SUBROUTINE SIC0CX,S,C> 

C SUBROUTINE TO EVALUATE SIN AND COS INTEGRALS 
IF(ABSCX) .GT. 1 .E-12)G0 TO 25 
S=0 .0 
00.0 
RETURN 

25 IF< ABS(X) *LT« 0 • 1 ) GO TO 26 
C SICI IS A SUBROUTINE IN CDC MATH-SCIENCE LIBRARY. 

CALL SICKS, C,X> 

S=S+1 .570 796 

C=. 57721 57+ALOGCABSCX) )-C 
RETURN 
S-X 
C=0.0 
N= 1 
NO 1 


26 
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NS- - 1 
NTC-0 
NTS=0 
FACT=1 . 

COUNT” 1 * 

PROD=X 
DO 10 1=2,10 
C0UNT=C0UNT+1 . 

FACT=FACT*COUNT 
PR0D=PR0D*X 
TERM=PROD/C COUNT* FACT) 
IFCN.EQ.DGO TO 15 
N=1 

IFCNS.EQ.-l > TERM* -TERM 

NS=-NS 

S»S+TERM 

I FCABSC TERM/S) *LT. 1 .E-5)NTS=1 
GO TO 112 
15 N=0 

IF(NC.EO.-t )TERM=-TERM 

NC=-NC 

C=C+TERM 

I F<ABS (TERM/C > .LT. 1 *E-5>NTC=1 
112 IFCNTS.EQ.O)GO TO 10 
IFCNTC.EQ. 1 ) RETURN 
10 CONTINUE 
RETURN 
END 

EOI ENCOUNTERED. 1 FILES. 

BYE 

ME650AC LOG OFF 10.37.33. 
ME650AC SRU 0.025 UNTS. 
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Table III-l. Summary of Rolling Experiments 
Film Thickness and Traction. 


V 

E 

V 

s 

P H 

t bath 

h 

c 

h 

m 

TC 

m/s 


m/s 

GN/m 2 

°c 

ym 

ym 



0.5 

2.00 

1.00 

0.81 

34 

.18 

.076 

.061 

ft 

0.60 

0.30 

ii 

37 

ff 

fl 

.056 

fl 

0.40 

0.20 

ff 

38 

VI 

IV 

.055 

IV 

0.20 

0.10 

If 

38 

If 

If 

.050 

11 

0.00 

0.00 

ii 

39 

ii 

ti 

.017 

11 

-.20 

-.10 

ff 

38 

ff 

ti 

— 

VI 

-.40 

-.20 

II 

40 

11 

If 

— 

ff 

-.88 

-.44 

VI 

41 

vv 

ff 

— 

0.75 

2.00 

1.50 

ii 

39 

.19 

.089 

.058 

m 

1.20 

0.90 

it 

ii 

.20 

fl 

.058 

ff 

0.67 

0.50 

ii 

it 

.20 

ti 

.055 

ff 

0.13 

0.10 

II 

ii 

.22 

ii 

.049 

n 

0.00 

0.00 

ii 

38 

.22 

ff 

.024 

it 

2.00 

1.50 

0.66 

33 

.22 

.15 

.051 

ft 

1.20 

0.90 

ff 

34 

.24 

.18 

.051 

ff 

0.67 

0.50 

ff 

34 

If 

ff 

.051 

II 

0.19 

0.14 

VI 

35 

If 

ti 

.051 

II 

0.00 

0.00 

ii 

VI 

ti 

ff 

.014 

If 

-.10 

-.08 

II 

ti 

•i 

•i 

.041 

it 

2.00 

1.50 

1.02 

40 

.15 

.064 

.061 

ff 

.27 

.20 

VI 

38 

.18 

.076 

.059 

fl 

0.00 

0.00 

If 

37 

.22 

it 

.020 

If 

-.27 

-.20 

vv 

ii 

.17 

ii 

.057 

0.51 

0.00 

0.00 

If 

If 

.18 

If 

.020 

0.72 

0.00 

0.00 

vv 

IT 

.22 

ti 

.023 

1.00 

0.00 

0.00 

ft 

If 

.24 

.127 

.023 

1.25 

0.00 

0.00 

11 

11 

.28 

.178 

.027 
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Table III-2. 


Summary of Rolling Experiments 
Ball Surface Temperature Rises, 
(P H =1.02 GN/m 2 ) 


V 

T. 

V 

s 

t bath 

CENTER 

AT 

MAX, 

AT 

m/s 


m/ s 

°c 

°c 

°c 

0.75 

2.00 

1.50 

42 

54 

64 

0.75 

1.20 

0.90 

42 

35 

36 

0.75 

0.99 

0.74 

43 

26 

26 

0.75 

0.56 

0.42 

42 

14 

15 

0.75 

0.27 

0.20 

40 

6 

6 

0.75 

0.13 

0.10 

41 

4 

4 

0.75 

0.00 

0.00 

41 

4 

4 

0.75 

-.27 

-.20 

40 

6 

6 

1.00 

2.00 

2.00 

44 

53 

62 

1.00 

1.50 

1.50 

45 

42 

46 

1.00 

1.30 

1.30 

45 

37 

39 

1.00 

1.04 

1.04 

44 

36 

37 

1.00 

0.60 

0.60 

43 

21 

21 

1.00 

0.40 

0.40 

43 

10 

10 

1.00 

1.12 

0.12 

42 

6 

6 

1.00 

0.00 

0.00 

41 

5 

5 

0.50 

0.00 

0.00 

41 

2 

2 

1.25 

0.00 

0.00 

41 

7 

7 
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Table IV-1. Frequency Shift and Velocity of Sound 
as Function of Pressure for 5P4E at 
24.4C (76F) (History B) . 


Date: December 26, 1975 

Separation Distance - 0.50 cm 

Resulting P = 24,700 psi 
8 






= 0.17 

GPa 

Spectrum 

Number 

Pressure 

kpsi* 

Frequency 

Shift 

-1 

cm 

Frequency 

Shift 

GHZ 

Refractive 

Index 

Velocity of 
Sound, m/s 

1 

4.78 

0.33194 

10.00 

1.6405 

2208 

2 

7.2 

0.34034 

10.21 

1.6453 

2258 

3 

10,460 

0.35364 

10.61 

1.6518 

2337 

4 

14,250 

0.36805 

11.04 

1.6594 

2421 

5 

16,800 

0.37819 

11.35 

1.6645 

2480 

6 

20,840 

0.39429 

11.83 

1.6726 

2573 

7 

23,125 

0.40259 

12.08 

1.6772 

2620 

8 

26,050 

0.41204 

12.36 

1.6829 

2672 

9 

28,200 

0.41642 

12.49 

1.6873 

2695 

10 

29,950 

0.42067 

12.62 

1.6909 

2715 

11 

31,500 

0.42347 

12.70 

1.6939 

2729 

* 

Multiply entry by 

6.894 X 10" 

•3 

to get pressure units of 

GPa. 
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Table IV-2. Frequency Shift and Velocity of Sound as a 

Function of Temperature at 0.40 GPa (59,000 psi) 
for 5P4E Oil (History A) . 


Date: March 2, 1976 

Separation Distance = 0.725 cm 
Resulting: T g = 70.3 C 


Spectrum 

Number 

Temperature 

C 

Frequency 

Shift 

-1 

cm 

Frequency 

Shift 

GHZ 

Refractive 

Index 

Velocity 
of Sound 
m/s 

64 

82 

0.43966 

13.19 

1.7222 

2787 

65 

80 

0.44236 

13.27 

1.7233 

2801 

66 

77 

0.44618 

13.39 

1.7246 

2823 

67 

75 

0.44973 

13.49 

1.7258 

2844 

68 

72 

0.65369 

13.61 

1.7271 

2867 

69 

70 

0.45511 

13.65 

1.7280 

2875 

70 

67 

0.45691 

13.71 

1.7292 

2883 

71 

65 

0.45789 

13.74 

1.7303 

2888 

72 

63 

0.45910 

13.77 

1.7313 

2894 

73 

59 

0.46064 

13.82 

1.7328 

2901 

74 

56 

0.46216 

13.86 

1.7343 

2908 
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Table IV-3. Least-Square Expressions for Sound Velocity, Temperatures and Sound 

Velocity at Transition for Different Constant Pressures: 5P4E History A. 


Formation 

Range of 


Least-Square 

Glass Transition 


Pressure, 

Measurements 

Expression for 

Temperature 

Sound Velocity 

MPa 

C 

Sound Velocity 

(T ) C 
g 

at T^, m/s 



v T 

_ 

2961 


7.99T 



210 

48-27 

L 





34 

2692 



V 

= 

2798 

_ 

3.15T 





g 









V T 

= 

3043 

— 

7.49T 



248 

49-27 

L 





41 

2739 



V 

= 

2856 


2.88T 





g 









V 7 

= 

3138 

- 

7.94T 



276 

58-32 

Li 





46 

2769 



V 


2903 


2.88T 





g 









V L 

= 

3172 

- 

6.79T 



321 

68-43 

Li 





53 

2812 



V 

3 

2942 

_ 

2.43T 





g 









V T 

= 

3259 

- 

6.95T 



355 

72-46 

L 





59 

. 2848 



V 

g 

3 

3015 

- 

2.48T 





V T 

— 

3410 

_ 

7.60T 



407 

82-54 

L 





70 

2875 



V 

3 

3032 

_ 

2.23T 





g 









Table IV-4. Experimental Fluids* In Volumetric 
Glass Transition Experiments. 


5P4E 
Nl, N2 
PI 

MCS 460 
MCS 1218 
Sanotrac 50 

Formulated Advanced Ester 


XRM-177-F4 
DN 600 
FN 2961 

52- Polybutene 
Sl-Diester 

53- DC 200-Blend I 
DC 200-50 


*Detailed descriptions are given in Appendix A. 
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Table IV-5. Fluids for Which No Transition was Observed 
by the Volumetric Method in the Indicated 
Temperature and Pressure Ranges (Isothermal 
Compression) . 


SAMPLE 

MAXIMUM PRESSURE 
(GPA) 

TEMPERATURE 

(°C) 

DIESTER (SI) 

1.13 

-3.3 

DC 200 BLEND II 

** .61 

7.2 

DC 200 - 50 (S3) 

* .92 

21.1 

ADVANCED ESTER 

*1.07 

16.9 


* PISTON SEIZED 
** DISPLACEMENT LIMITED 
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Table IV-6. Approximate Yield Shear Strength 
at Atmospheric Pressure 


FLUID 

T 

T 

g 

- t yield, KPa 

5P4E 

- 60C 

- 23C 

5.38 

N1 

< T 

g 

- 62C 

3.24 

XRM 

< T 

g 

N.A. 

2.96 
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IR MICRODETECTOR OR 
OPTICAL MICROSCOPE 
AND CAMERA 


SAPPHIRE 


ROTATING BALL 

COUNTER- 
BALANCE ^ 


SAPPHIRE SUPPORT 
ASSEMBLY 


-OVERFLOW PAN 

fl 


SUPPORT PADS 
(SAPPHIRE) 
CONTACT FEED 
RESERVOIR 

NORMAL FORCE LOAD 
CELL 

SUPPORT BEAM 


HEATING COIL 


CONSTANT TEMPERATURE BATH 



VARIABLE 
SPEED MOTOR 


— CONTACT 
LOADING 
: WEIGHTS 

AIR BEARING 


OIL PUMP MAIN RESERVOIR 


Figure III -I. Schematic of Sliding EHD Contact Simulator. 
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GEARING 
VARIABLE SPEED 
D.C. MOTOR 


OVERFLOW PAN 


CONSTANT 
TEMPERATURE 
BATH 


IR MICRODETECTOR 
or 

OPTICAL MICROSCOPE 
AND CAMERA 



THRUST BEARING 
ROTATING SAPPHIRE 
LOADING WEIGHTS 

SUPPORT BEAM 
SUPPORT COLUMN 


VARIABLE SPEED 
D.C. MOTOR 


OIL PUMP 


OIL RESERVOIR 


Figure III-2. Schematic of the Combined Rolling and Sliding EHD Contact Simulator. 







BALL SURFACE TEMPERATURE RISE AT 


NON-DIMENSIONAL PARAMETER L FOR P u = 1.70 GPa. 

H 

.1 .5 1 5 10 50 100 



Figure 1 1 1-3. Ball Surface Temperature Rise at Contact Center (.011 pm Ra roughness, 
1 < A < 2 for P H = 2.03 GPa, otherwise A > 2). 




Figure 1 1 1-4. 




Surface Profiles in Direction of Sliding 
(P H = 1.24 GN/m^, V $ = 1.02 m/s. Hertz 
Diameter d = 0.43 mm, pure sliding). 




oo 2d d . 5d .3d .2d .15d . 1 2d .10d .08d X 


a) Unused ball surface. 



- 2d d . 5d .3d .2d .15d . 1 2d .10d .08d X 

b) Ball surface after running 24 minutes. 


Figure 1 1 1-5. Power Spectra of the Surface Profiles 
(Pn = 1.24 GN/mS V = 1.02 m/s. Hertz 
diameter d = 0.43 mm, pure sliding.) 
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0 po^ition 2 ^hift 3< ? b 4d 5d 


a) Unused ball surface (.38 ym Ra roughness). 



POSITION SHIFT , 3 

b) Ball surface after running 24 minutes. 


Figure III-6. Autocorrelation Functions of the Surface Profiles. 

(.X > 3d not included, P„ = 1.24 GN/nr, Vc = 1.02 m/s. 
Hertz diameter d = 0.43 mm, pure sliding) 
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Radiation 
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.2 .4 .6 

Fluctuations after Running 8 Minutes, 


...... 



Figure III-7. Normalized Histograms of Ball Surface Radi atipn Fluctuations 
(.38 ym Ra initial roughness, P 
Hertz diameter d = 0.43 mm, pur 


H = 1.24 GN/nr , V c = 1.02 m/s, 
e sliding). 






b) Ball Surface Temperature Fluctuations after 8 Minutes Run-in. 

Figure III -8 . Power Spectra of Ball Surface Temperature Fluctuations 
(.38 pm Ra initial roughness, P., =1.24 GN/m , 

V = 1.02 m/s. Hertz diameter -0.43 mm, pure sliding). 
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POSITION ON FERROGRAM 





FILM THICKNESS h, m 




Figure III-10. Film Thickness as a Function of Slide-Roll 
Ratio (.011 pm Ra roughness). 
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TRACTION COEFFICIENT TC 



- 1.0 0 0 .5 1.0 1.5 2.0 

i 'OF -ROLL RATIO >; 

Flqurf’ ' ■ ! ' i i'i.i Coefficient as a Function of Slide-Roll 

R - : ■ 1 ■ • (.nil uni Ra roughness). 
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FILM THICKNESS h, ym 



Figure III-12. Film Thickness and Traction Coefficient as a Function of Rolling 
Velocity (£ = 0, .011 urn Ra Roughness, P^ = 1.02 GPa). 


TRACTION COEFFICIENT TC 



TEMPERATURE RISE 



SLIDE-ROLL RATIO I. 


Figure 1 1 1 - 1 3 . Ball Surface Temperature Rise as a Function of 

Slide-Roll Ratio (.011 urn Ra Roughness, P M = 1.02 GPa , 
V = 0.75 m/s). H 
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TEMPERATURE RISE AT 




TEMPERATURE RISE AT 



ROLLING VELOCITY 9, m/s 

Figure III-15. Ball Surface Temperature Rise as a Function of Rolling Velocity (.011 pm Ra 
Roughness, P H = 1.02 GPa, E = 0, Maximum AT^ Occurs at Contact Center). 



TEMPERATURE RISE 





TEMPERATURE RISE A 



DISTANCE FROM CENTER, mm 


Figure HI-17. Ball Surface Temperature Rise Along the Contact 
Centerline (.011 urn Ra Roughness, P u = 1.02 GPa, 
V = 1.00 m/s) . H 
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-.25 -.20 -.15 -.10 -.05 0 .05 .10 .15 .20 .25 

INLET EXIT 

DISTANCE FROM CENTER, mm 


Figure III-18. Ball Surface Temperature Rise Along Contact 

Centerline (.011 wm Ra Roughness, P u = 1.02 GPa, 

>; = 0). H 
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PREDICTED AVERAGE TEMPERATURE RISE, 



Figure 1 11-19. Comparison of Predicted Average and Actual 
Maximum Ball Surface Temperature Rises for 
L > 5. Combined sliding and rolling as pe 
Equation 21. 
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Figure I V-l . Thermal Expansion Coetr 1 - v- Heat 
for a Typical Glassy R-; i t : •, 1 . 
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SPECIFIC VOLUME 



Figure IV -2 - 


Typical Plot of Specific Volume Versus 
Temperature for a Glassy Material at 
Different Formation Pressures. 
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» 


SPECIFIC VOLUME 





Tgs Tg 2 Tg, 

temperature > 


Figure IV-3. Dependence of the Glass Transition Temperature 

on the Cooling Rate, . L = Liquid and G = Glass. 
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Figure IV-5. Frequency Spectrum of a Naphthenic Base Oil 
at 24. 4C and Atmospheric Pressure. 



Figure IV-6. The Variation of the Brillouin Frequency Shift and the Sound Velocity 
with Temperature at Atmospheric Pressure for Atactic Polystyrene [53] 


Sound Velocity-m/s 






Figure I V-7. Schematic Arrangement of the 
Light-Scattering Experiment. 
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i 


Corner 



Figure IV-8. Schematic Arrangement of the Fabry-Perot Interferometer. 



Transmission 
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Windows 


Figure IV-10. 


r^n The Intensifies 
•fmachef PerSndicuUri^the Paper Plane at the 
Centerline of the Cel . 
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Figure IV-11 . Schematic of Formation Histories Used to Form the Glass 
A. and B Represent Constant Formation Histories at 
Pressure P.. and Room Temperature Respectively. 
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Figure IV-12. 


Frequency Spectrum for 5P4E Fluid Recorded at 0.19 GPa and 24. 4C (History B) 
B and R represent the Brillouin and Rayleigh Components Respectively. 



Figure IV-13. Variation of Frequency Shift and Velocity of Sound 
with Pressure at 24. 4C (76 F) for 5P4E Fluid 
(History B). Arrows Indicate Glass Transition. 
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Frequency Shift 



Temperature , c 


Figure IV-15. Variation of Frequency Shift and Velocity of Sound 
with Temperature at 0.40 GPa (59,000 psi) for 5P4E 
(History A). Arrows Indicate Glass Transition. 
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i 


Velocity Of Sound , m/s 



Frequency Shift (« 








99T 

temperature- 



Figure IV-18. Phase Diagram for 5P4E Fluid Based on Light Scattering Experiment. = 183 C/GPa. 


History A. History B. 




Figure IV-19. High Pressure Dilatometer (.13 to 1.2 GPa). 
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Figure I V -20 . Low Pressure Dilatometer. 
(.014 to .27 GPa) 


166 






167 


OnH 
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P (GPa) 


Figure IV-26. Dispersion Region in Isothermal Compression by Dilatometry. 
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LOG-, n T(MPa) 



Figure IV-29. Shear Stress Between Concentric Cylinders for 5P4E 
Cooled at 2.75 C/min., Shear Rate ? = .435 sec" 1 . 
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5P4E 



SLIDING 






APPARENT VISCOSITY, N-sec/m' 



PRESSURE (GPa) 


The Variation of Apparent Viscosity with 
Pressure and Temperature. 



811 

TEMPERATURE (°C) 



<$> REFERENCE [3] DATA 

<>♦ REFERENCE [2] DATA (ONSET AND 
- MIDPOINT RESPECTIVELY) 

LIGHT SCATTERING HISTORY B 
D LIGHT SCATTERING HISTORY A 
O PVT HISTORY A 
• PVT HISTORY B 



RELATIVE THRUST 






BALL SURFACE TEMPERATURE 



Figure IV-34. The Physical State of N1 Lubricant in the Center 
of a Sliding Contact at Various Speeds and Hertz 

Pressures (Reference [9]) and 

Represent the Light-Scattering 

and the Volumetric Glass Transition Lines 
Respectively. 
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Ball Surface Temperature, 



Sliding Speed, m/s 


Figure IV-35. The Physical State of N1 Lubricant in the Center of a Sliding Contact 
at Various Speeds and Hertz Pressures and — Represent the 

Light-Scattering and the Volumetric Glass Transition Lines. 




(a) Light Scattering 
Transition Data 



1 = 2.0 
1 = 1.2 
1=0.99 


(b) Volumetric 
Transition Data 



1 = 2.0 
Z= 0.99 
1=0.56 
1 = 0.27 


Figure IV-36. Ball Surface Temperature Distribution Along the Contact 
Centerline at Constant Hertzian Pressure of 1.0 GPa 
and Constant Rolling Speed of 0.75 m/s. (a) and (b) 
Represent the Physical State of the Lubricant in the 
Contact Area as Determined by the Light Scattering and 
Volumetric Transition Data Respectively. 
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X 


a) Symmetric Pressure Distribution. 




b) Antisymmetric Pressure Distribution. 


Figure V-3. Decomposition of the EHD Pressure-Cylinder. 
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x 



a) Symmetric Pressure Distribution. 


y 



b) Antisymmetric Pressure Distribution. 


Figure V-4. Decomposition of the EHD-Pressure-The Plate. 


185 





Figure V-5. Schematic Diagram of EHD Contact Zone. 
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Figure V- 6. Test Case 1 - Uniform Pressure Distribution. 



P hz : Maximum Hertz Pressure Hertzian 



Figure V-7. Test Case 2 - Hertzian Contact. 


188 



Figure V-8. Test Case 3 - Polynomial Pressure Distribution. 
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Figure V-9. Test Case 3 - Random Errors Introduced into Film Thickness Values. 
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Figure V-10. Test Case 3 - Errors Introduced in Location of Inlet Point. 
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Figure V - 1 1 . Calculations Based on Experimental Film 
Thickness Measurements [71]. 


192 


I 


APPENDIX A 


DESCRIPTION OF EXPERIMENTAL FLUIDS 


The following table summarizes the oils investigated in this 
research and gives characteristic data for each oil. 

Experimental Fluids 

Description 
Polyphenyl Ether 
Cycloaliphatic Hydrocarbon 
Synthetic Hydrocarbon 


Symbol 

5P4E 

MCS-1218 
MCS-460 
Santotrac 50 


Cycloaliphatic Hydrocarbon Traction Fluid 
plus additive 


N1 


Naphthenic Base Oil R-620-15 


PI 

Nap 2 
FN 2961 
XRM-177-F4 

DN-600 

SI 

DC-200-Blend II 
S3-DC-200-50 


Paraffinic Base Oil R-620-12 

Napthenlc Base Oil R-620-16 

Super Refined Napthenlc Mineral Oil 

Synthetic Paraffinic Hydrocarbon 
plus anti-wear additive 

Polyalkyl Aromatic 

Diester 

Dimethylsiloxane Blend 
Dimethylsiloxane 


Advanced Ester 


S2 


Polybutene LF-5193 
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Fluid Characterization 


Symbol : 
Type: 

Source : 
Properties : 


5P4E 

Five-ring Polyphenyl Ether 
Monsanto Company 

2 

Viscosity at 37. 8C, m / s 
Viscosity at 98. 9C m 2 /s 
Density at 22. 2C, Kg/m 3 
Density at 37. 8C, Kg/m 3 
Flash Point, C 
Pour Point, C 


363 x 10 
13.1 x 10" 6 
1.205 x 10 3 
1.19 x 10 3 
288 
4.4 


Symbol : 
Source : 
Type: 

Properties : 


MCS-1218 

Monsanto Company 
Cycloaliphatic Hydrocarbon 

It is a combination of two components each have a 
molecular weight less than 1000. 

2 

Viscosity at 37. 8C, m /s 1418 x 10 

Viscosity at 98. 9C, m 2 /s 18.37 x 10 

Density at 23. 9C, Kg/m 3 0.94 x 10 3 
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Symbol: 

Source: 

Type: 

Properties : 


MCS-460 

Monsanto Company 

Synthetic Hydrocarbon 

2 

Viscosity at 37. 8C m /s 

Viscosity at 98. 9C m 2 /s 

2 

Viscosity at 148. 9C m /s 

Pour Point 

Density 25C Kg/m 2 


37.2 x 10 
4.0 x 10~ 6 
1.9 x 10 -6 
-29 to -32C 
.9327 x 10 3 


Symbol : 
Source : 
Type: 

Properties : 


Santotrac 50 

Monsanto Company 

Synthetic Cycloaliphatic 

2 

Viscosity at 37. 8C m / s 

2 

Viscosity at 98. 9C m /s 
Pour Point C 
Density at 37 . 8C 
Flash Point C 
Fire Point C 


Hydrocarbon Traction Fluid 
.34 x 10 -4 
.056 x 10~ 4 
-37 

.889 x 10 3 Kg/m 3 

163 
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Specific Heat at 37. 8C J/Kg-K 0.51 

Additive package includes: Antiwear (zinc dialkyl 
dithiophosphate) , Oxidation inhibitor Antifoam, VI 
Improver (Polymethacrylate) . 
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Symbol : 

N1 


Source : 

Sun Oil Company 


Type: 

Naphthenic Base Oil R-620-15 


Properties: 

2 

Viscosity at 37. 8C, m /s 

24 x 10~ 6 


2 

Viscosity at 98. 9C, m / s 

3.728 x 10~ 6 


Viscosity Index (ASTM D-2270) 

-13 


Flash Point, C 

157 


Pour Point 

-43 


Density at 20C, gm/cc Kg/m 3 

0.9157 x 10 3 


Average Molecular Weight 

305 

Symbol : 

PI 


Source: 

Sun Oil Company 


Type: 

Paraffinic Base Oil 


Properties: 

Viscosity at 37. 8C m 2 /s 

33.74 x 10~ 6 


Viscosity at 98. 9C m 2 /s 

5.402 x 10‘ 6 


Density at 25C Kg/m 3 

.8602 x 10 3 


V.I. (ASTM D-2270) 

103 


Pour Point C 

-15 


Average molecular weight 

404 
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Symbol: 

Nap 2 


Source : 

Sun Oil Company 


Type: 

Napthenic Base Oil R-620-16 


Properties: 

? 

Viscosity at 37. 8C m /s 

114.2 x 10" 6 


2 

Viscosity at 98.9C m /s 

8.076 x 10” 6 


Viscosity Index (D2270) 

< 0 


Density at 20C Kg/m'* 

.9303 x 10 3 


Average Molecular Weight 

357 


Refractive Index 

1.5173 


Pour Point C 

-23 

Symbol : 

FN 2961 


Source: 

Humble Oil and Refining Company 


Type: 

Super Refined Napthenic Mineral Oil 


Properties : 

Viscosity at 37. 8C m 2 /s 

78.08 x 10~ 6 


Viscosity at 98. 9C m 2 /s 

8.238 x 10* 6 


2 

Viscosity at 148. 9C m /s 

3.3 x 10“ 6 


Density at 15. 6C Kg/m 3 

.887 x 10 3 


Pour Point C 

-34 
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Symbol: 


XRM-177-F4 


Source: Mobil 

Type: Synthetic Paraffinic Hydrocarbon 

plus antiwear additive 

2 

Properties: Viscosity at 37 . 8C Ns/m 

Viscosity at 98. 9C Ns/m^ 

Pour Point C 
Density 37. 8C Kg/m 3 


376 x 10 J 
31.6 x lCf 3 
< -40 

.8389 x 10“ 3 


Symbol: DN-600 

Source: Continental Oil Company 

Type: Polyalkyl Aromatic 

2 

Properties: Viscosity at 37. 8C m /s 

2 

Viscosity at 98. 9C m /s 

2 

Viscosity at 148. 9C m /s 
Density at 37. 8C Kg/m 3 
Pour Point C 

Specific Heat at 37. 3C J/Kg C 


30 x 10 
5.0 x 10* 6 
2.3 x 10~ 6 
0.851 x 10 3 
- 60 
1624 
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Symbol : 
Source: 
Type: 

Properties : 


Symbol : 
Source : 
Type: 

Properties : 


Sl-Dies ter 

Rohm and Haas Company 

Diester-Plexol 201 bis-2-ethyl hexyl 

2 

Viscosity at -53. 9C m /s 

? 

Viscosity at 37. 8C m / s 

2 

Viscosity at 98. 9C m /s 
Cloud Point (ASTM D-2500) 


sebecate (PL 5159) 
7988 x 10“ 6 
12.75 x 10“ 6 
3.32 x 10~ 6 
below -54C 


DC-200-Blend II 

Dow Corning Corporation 

Dimethylsiloxane Blend E1923-49 


Viscosity at 25C m /s 1060 x 10 

Molecular Weight (GPC) 160,000 

Component Viscosities: 

50 x 10' 6 m 2 /s 75% 

11.000 x 10" 6 m 2 /s 12% 

330.000 x 10~ 6 m 2 /s 13% 
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Symbol: 

Source: 

Type: 

Properties : 


S3-DC-200-Blend II 

Dow Corning Corporation 

Dimethyls ilixane DC-200 

? 

Viscosity at 37. 8C m /s 

2 

Viscosity at 98. 9C ra /s 

3 

Density at 25C Kg/m 
Molecular Weight (Average) 


82.6 x 10 
33.1 x 10“ 
.968 x 10 3 * 
7,000 


Symbol: 
Source : 
Type: 

Properties : 


Advanced Ester 

Shell Oil Company 

Based on Pentaerythritol 
Base Oil) 

2 

Viscosity at 37. 8C m /s 

2 

Viscosity at 98. 9C m /s 

2 

Viscosity at 148. 9C m /s 

3 

Density at 37. 8C Kg/m 


(Aeroshell Turbine Oil 555 

25.8 x 10" 6 
5.1 x 10~ 6 
2.3 x 10“ 6 
.979 x 10 3 


200 


Symbol: 

Source: 

Type: 

Properties: 


S2 

American Oil Company 
Polybutene LF-5193 

? 

Viscosity at -17. 8C m /s 

2 

Viscosity at 37. 8C m / s 

Viscosity at 98. 9C m 2 /s 

Density at 25C Kg/m 3 

V.I. (ASTM D-2270) 

Polymer Number Average molecular 
weight 


18,836 x 10 
109 X 10~ 6 
10.6 x 10 ^ 
0.8443 x 10 3 
87 

409 
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Elastrohydrodynamic Forces", M.S. Thesis, Georgia Institute of 
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